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The statistical mechanical treatment previously applied 
to homogeneous polymer solutions has been extended to 
heterogeneous polymers composed of numerous molecular 
species differing in polymer chain length. Entropies of 
mixing with solvent are derived with reference to three 
standard states: the oriented pure components, the dis- 
oriented pure components, and the disoriented mixture of 
species. The entropy of mixing with respect to the last- 
mentioned state is identical with that previously derived 
neglecting polymer heterogeneity. Upon introducing a 
van Laar heat of mixing term, expressions for the partial 
molal free energies are derived. These should be reasonably 
correct in view of the satisfactory agreement between 
observed solvent activities in polymer-solvent solutions 
and theory, as shown-by Huggins. The relationships should 
be useful not only in dealing with solutions, but also in 
equilibria involving heterogeneous polymers in the absence 
of solvent. The thermodynamic equations are applied to 
the problem of equilibrium in the range of partial misci- 


bility for solvent-heterogeneous polymer systems. The 
requisites for efficient fractionation of high polymers are 
discussed. The efficiency of separation depends on the 
ratio of the volumes of the supernatant and the precipitated 
phases. In order to attain a high value for this ratio, very 
dilute solutions must be employed. The higher the molecu- 
lar weight at which separation occurs the greater the 
dilution required for the same sharpness of separation. 
The free energy expressions are also applied to the formu- 
lation of equilibrium constants to be employed in reversible 
polymerization-degradation processes. It is shown that 
concentrations should be expressed in moles per unit 
volume in the equilibrium constant. The equilibrium 
state for a polyfunctional condensation, considered as a 
reaction between functional groups, is not affected by the 
polymeric nature of the reacting species; the position of 
equilibrium should be the same as is found for analogous 
monofunctional reactants under the same conditions. 





INTRODUCTION 


OLUTIONS of high polymers have been 
known for some time to behave anoma- 
lously?* with respect to the ideal solution laws 
as originally set up for solutions of simple mole- 
cules. These apparent anomalies in thermo- 
dynamic and related behavior have been largely 


cleared up by statistical mechanical treat- 


1 A preliminary report on this work was published pre- 
viously, J. Chem. Phys. 12, 114 (1944). 

2 W. Ostwald, Kolloid Zeits. 46, 248 (1928). 

’K. H. Meyer and R. Liithdemann, Helv. Chim. Acta, 
18, 307 (1935). 

*K. H. Meyer, Natural and Synthetic High Polymers 
iapatanne Publishers, Inc., New York, 1942), pp. 582- 


ments®~ of high polymer solutions which have 
appeared in the last few years. These treatments 
are based on the following physical approxima- 
tion: The liquid mixture composed of long chain 
polymer and monomeric solvent molecules is 
assumed to conform to a pseudo-lattice structure, 
each cell of which may be occupied either by a 
molecule of solvent or by an element of the 
chain (segment) equal in size to a solvent mole- 
cule. Proceeding from this concept, it has been 


nN 


5M. L. Huggins, J. Phys. Chem. 46, 151 (1942); Ann. 
to ieee Sci. 43, 1 (1942); J. Am. Chem. Soc. 64, 1712 
‘P. J. Flory, J. Chem. Phys. 10, 51 (1942). 

7A. R. Miller, Proc. Camb. Phil. Soc. 39, 54 (1943). 
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possible to derive by an approximate mathe- 
matical procedure the following expression® for 
the entropy of mixing ” solvent with N polymeric 
solute molecules (initially in random configura- 
tions in the pure state) where 


AS mixing = — RE In v1: +N In v2]; (1) 


v; and v2 are the volume fractions of the com- 
ponents, solvent and polymeric solute, respec- 
tively. This expression differs from the one 
representing the so-called ‘“‘ideal’’ solution in 
that mole fractions are replaced by volume frac- 
tions. Inasmuch as the ideal solution law can 
be derived theoretically only for solutions of 
equal-sized species, where volume fractions equal 
mole fractions, the conventional ideal solution 
is seen to be included as a special case of Eq. (1). 

Although certain limitations on this statistical 
treatment as applied to polymer solutions have 
become apparent,® sufficient agreement with 
experiment has been obtained® ®*!° to encourage 
extension of this method of attack. It is the 
purpose of the present paper to extend the 
previous statistical mechanical method to hetero- 
geneous polymers and their solutions in which the 
polymer constituent is actually a mixture of 
many species differing in length, or molecular 
weight. The resulting equations afford a basis for 
dealing with chemical equilibria involving inter- 
conversion of the polymer species. They can also 
be applied to the problem of fractionation of 
polymers, as, for example, by fractional pre- 
cipitation from solution. 


DERIVATION OF ENTROPIES OF MIXING 


The procedure closely parallels that which has 
been given previously for a polymer consisting of 
a single species. The pseudo-lattice for the solu- 
tion is to be filled by successively adding the 
polymer molecules to it, then filling in the vacant 
cells with molecules of solvent. Suppose that j —1 
polymer molcules previously have been assigned 
locations in the lattice. The number of unoc- 


8 Paper presented before the Division of Physical and 
Inorganic Chemistry at the Pittsburgh Meeting of the 
American Chemical Society, September 6, 1943. To be 
published shortly. 

9M. L. Huggins, Ind. Eng. Chem. 35, 216, 980 (1943); 
Ann. N. Y. Acad. Sci. 44, 431 (1943). 

10 G. Gee and L. R. G. Treloar, Trans. Faraday Soc. 38, 
147 (1942). 
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cupied sites in the lattice will be given by 
j-1 
No— Do Xi 
n=1 


where mp is the total number of lattice sites and 
x; is the number of segments in the ith molecule 
(t<j). This factor will represent the number of 
sites available to the terminal segment of the 
jth molecule to be added to the lattice. The 
cell to which this first segment is assigned will be 
surrounded by y-sites for the second segment, 
where y is a “coordination number’’ for the 
lattice. Some of these may be occupied by seg- 
ments of molecules previously added to the 
lattice. The expected number of available sites 
for the second segment of the jth molecule can 
be expressed therefore as the product of y into 
the fraction of all lattice sites which are not 
occupied by previously added segments of mole- 
cules ;* this product is given by 


j-1 
y{ nm—->d v1) [me 
i=l 


Of the y-lattice sites surrounding the second 
segment, one is occupied by the first segment. 
The expected number of sites available to the 
third segment, therefore, will be 


(y—1) mo—E 21-2) /'me 


Similarly for the th segment of the jth molecule 
the expected number of available sites will be 


(y—1)| mE v0] /m (2)" 


The number (or expected number) of con- 
figurations available to the jth molecule as a 
whole can be obtained by combining the factors 
given above for the expected numbers of sites 
available to each segment. The resulting ex- 


1 The term (—1) in these expressions may be open to 
objections on the grounds that segments of the jth molecule 
are not to be considered the equivalent of segments of other 
molecules which may locate any place in the lattice. In any 
event, its contribution to the term is insignificant, and 
could be omitted entirely. It has been included here for the 
reason that it simplifies reduction of the subsequent 
equations. 
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pression becomes 


1 j-1 
ra me x «:|E4/tr- 1)] 


xI | (m-E m—e41 or /no} 


The symmetry number o of the polymer mole- 
cules usually may be taken equal to two, in con- 
sideration of the equivalence of the ends of the 
molecule. Replacing y/(y—1) by unity in the 
second factor in brackets in the equation given 
above, 


zj—1 


vj=[mo/o(y—1) ] Tt 
of (-E-1)r—n fn] 
“(F) (»-E=)'/ 
(10 z x) 1, (3) 


The total number of configurations available 
to the system of N polymer molecules is given by 


w= ( 1/TL Ne! IT > (4) 


j=1 


where N, is the number of polymer molecules 
composed of x units. Substituting Eq. (3) in (4) 
and simplifying the product term in the nu- 


merator, 
N 


w-[1/(o tN Dy 
| net / m—-¥ x)}] 


N N @ 
LD xj=L x1=L xN,, 
im] 


j=1 z=1 


Noting that 


and replacing mo with 


no=n+)>. xN, 


z=1 





where x is the number of solvent molecules, 
W=[1/(o* TT Nz!) L(y¥—1)/(n+ Dx.) Pe D9 
X[(n+>¥ xN,)!/n!]}. (5) 


In Eq. (5) and in those which follow, the limits 
on the summations are omitted inasmuch as they 
invariably run from x=1 to ~. 

The entropy of mixing is given by k ln W, 
where & is Boltzmann’s constant. Or, if at this 
point » and the N,’s are expressed in moles, 
rather than as numbers of molecules, 


AS mizing = R In W. 


Substituting for W from Eq. (5) and introducing 
Stirling’s approximation for factorial terms, 


n 
AS mixing = —Rjn In (——) 
7 n+ DxN, 


+z[W. ™ (=) 


~FL(x-N- En (y-1)-1}+W In o| 6) 
or 
ASmixing= —R{n In 11+ CN, In (v,/x) J 
—>[(x-1)N.][In (y—1)—-1]+NIno}, (6’) 


where v7; and v, are volume fractions of solvent 
and x-mer, respectively. In writing these equa- 
tions the term for solvent (first term) could have 
been included with the first summation, the 
solvent being considered as a ‘‘polymer’’ com- 
posed of one segment. However, in consideration 
of possible heats of interaction between solvent 
molecules on the one hand and polymer segments 
on the other, it will facilitate formulation of free 
energy expressions to keep these terms separate. 

In deriving Eqs. (6) and (6’) the liquid 
structure has been approximated by a lattice 
composed of cells equal in size to one solvent 
molecule. This physical approximation has re- 
ceived due consideration in previous publica- 
tions. In the case of heterogeneous polymers, 
however, it leads to a summation term” 
> N.1n (vz/x) which contains only integral 
values of x. In applying these and subsequent 

The other summation terms in Eqs. (6) and (6’) 


gen merely total number of segments and of mole- 
cules. 
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equations to actual heterogeneous polymers it 
will be necessary to let this summation include 
the actual molecular species, for which x, the 
ratio of the molecular volumes of polymer mole- 
cule and solvent, will in general be non-integral. 
The additionai approximation introduced by this 
extension will be negligible so long as the x values 
concerned are large. It should be borne in mind, 
however, that in subsequent equations the x 
values may not necessarily be integral. 

The above expressions (6) and (6’) refer to 
the entropy change on forming the solution from 
Ni, No, N3, etc. molecules of monomer, dimer, 
trimer, etc., each initially in an ordered (crystal- 
line) state, and nu molecules of solvent. This 
entropy change can be separated into three com- 
ponents: (a) the entropy of disorientation of the 
individual pure components, (b) the entropy of 
mixing these, and (c) the entropy of dilution of 
the mixture. 


AS nisies - ASa+tAS+AS.. (7) 
The quantity AS, is, of course, the sum of the 


entropies of disorientation of the individual 
polymer species.® For the x-mer species 
ASa, z= RN,[In (x/o) +(x —1) In (y—1) — (x—1)]. 
Hence, 
ASg=R [Nz In (x/o) ] 
+R YL (x—-1)N2JLIn (y—-1)—-1]. (8) 
Subtracting this entropy of disorientation 
from Eq. (6) and setting »=0, 
AS,=—-R>[N. In (xN./>,xN;z)]. (9) 
This equation represents the entropy of mixing 
together the previously disoriented individual 
polymer species. 
Finally, for the entropy change on diluting 


the mixed, disoriented polymer molecules with 
n molecules of solvent 


AS. _ AS ntctes — AS, = AS, 


= -R|n In baton) 
+2[v.6 (Felt « 
= —R(n|1nv,+N Inv), (10’) 
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where v2 is the volume fraction of all polymer 
molecules; i.e., 


Vo= dz. 


Equation (10’) is identical with the equation 
previously derived® for the entropy of mixing N 
disoriented polymer molecules, all of the same 
size, with 2 solvent molecules. When dealing with 
situations in which no change in the distribution 
of polymer species is encountered, this equation 
may be employed as the basic entropy of mixing 
expression without reservations for heterogeneity. 
In the expression previously derived® ¢ for partial 
molal quantities it is only necessary to replace 
the number of segments x with its number aver- 
age £,. In short, heterogeneity can be disregarded, 
provided that the polymer is confined to one 
phase and that no polymerization and/or de- 
gradative steps are involved within that phase. 
For example, vapor pressures and osmotic pres- 
sures of polymer-solvent mixtures can be handled 
using AS., as previously published,®®* without 
regard for heterogeneity of the polymer. It is 
necessary, however, that the number average 
molecular weight, rather than some other aver- 
age, be employed. 

When two phases containing polymer are in 
equilibrium, passage of polymer from one phase 
to another being permissible, the undiluted, 
heterogeneous mixture of polymer molecules of 
given distribution is no longer adequate as a 
reference state. For the treatment of certain 
problems, e.g., chemical equilibria between 
polymer species (i.e., condensation or addition 
reactions between polymer species, or equilibrium 
degradation of the polymer molecules), Eqs. (6) 
or (6’) should be employed as the basic entropy 
of mixing expressions. The state in which the 
individual species are separate and oriented then 
represents the standard reference state. 

In other cases it may be more convenient to 
employ the unmixed, disoriented species for the 
reference state. This reference state is sufficient 
for the treatment of the equilibria encountered 
in the fractionation of high polymers by frac- 
tional precipitation or solution. The entropy of 
mixing the separate, disoriented species with 
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solvent is given by 


an = AS,+AS, 


n 
= =R} nin (——} 
n+>-xN, 
xN, 
N,z In | ————_ 
+z| , bee we 
= —R[nInv,+D(N, Inv.) ]. (11’) 


An equation of this form for the entropy of 
mixing in a heterogeneous polymer solution was 
predicted by Gee!* as an empirical extension of 
the entropy of mixing expression® for a solution 
composed of a single polymer component. 


PARTIAL MOLAL ENTROPIES 


As has been pointed out above, depending on 
the application, any one of three standard refer- 
ence states—the disoriented mixture of polymer 
species in the absence of solvent ; the disoriented, 
separated polymer species; or the oriented pure 
components—may be preferable. If the first- 
mentioned reference state is to be employed 
(AS.) the partial molal functions previously 
derived,®>® wherein x is replaced by Zn, will 
suffice. For the other two cases either AS* of 
Eq. (11) or AS of Eq. (6) must be differentiated 
to obtain the partial molal entropies. These are 
given below. 

For the partial molal entropy of the solvent 
per mole, 


AS, = dAS/dn=AS,* = dAS*/dn 





=—R{In (—-—)+1 
n+> xN, 
| n > N-z 
n+>xN, n+ >Dx«N, 


= —R[In (1-2) +02—D(v2/x)]. 


Or, since the number average molecular size Z, 
is given by 


f,=> xN,/D N,=02/¥(02/x) 


AS,=A8,*= —R{In (1 —v2)+v2(1 _ 1/Zn) |, 


8G. Gee, Rubber Chem. Tech. 17, 672 (1944). 


which is identical with the expression previously 
derived,® * except that x has been replaced by Z,. 
For the polymer constituents 


AS,=dAS/dN,; 


N;z 
= -R{In (——_) #1 
n+>> xN,z 


xn x> N,; 
n+>xN, n+>xN, 





—(x—1)[In (y—1)—1]+1n “| 
= — R[In (v,/x) +02x(1—1/2,) 
—(x—1) In (y—1)+Ino]. (13) 
Similarly 
AS,* = —R[In v,—(x—1)+0ox(1—1/2n)]. (14) 


PARTIAL MOLAL FREE ENERGIES 


Assuming that the heat of mixing of solvents 
with polymer can be expressed in the van Laar 
form, then 


AH=B(Vyn) > V.N./(Vint+d V-N:z), 


where V; is the molar volume of solvent, V, is 
the molar volume of the x-mer polymer species, 
B is a constant characteristic for each solvent- 
polymer pair, and the summations are over all 
polymer species. This can be written 


AH=BV4n > xN./(n+>d xN-) (15) 


since V,=xVj. 
The partial molal heats of mixing are obtained 
by differentiation. 


AH, =0AH/dn=BV w?, 
AH,=0AH/dN,=Bx Vw. 


(16a) 
(16b) 


Combining these experiments with those given 
above for the partial molal entropies, 


AF, =AF,*=AH,—TAS, 
=RT[In (1—v2)+02(1—-1/%,) +027], (17) 
where » = BV,/RT. Similarly, 
AF,=RT[In (vz/x) +02x(1—1/Zn) 
—(x—1) In (y—1)+1n o+px(1—v2)?] (18) 
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and 
AF ,* =RT[In 9,—(x—1) +00x(1—1/n) 
+pyx(1—ve)?]. (19) 


Huggins® *® has shown that the data available 
on solvent activities in solvent-polymer mix- 
tures agree very well with Eq. (17) when yu is 
suitably chosen. In accordance with theory, u 
is found to be characteristic for a given solvent- 
polymer pair at a given temperature, quite 
independent of concentration. Huggins has 
tabulated » values for various systems. 

Only in one case have solvent activities been 
measured carefully over the entire concentration 
range, namely, the benzene-rubber system thor- 
oughly investigated by Gee and Treloar.!° 
Huggins has shown that a single value of yu suf- 
fices to establish a good correlation between Eq. 
(17) and observed activities at all concentrations. 
Gee and Treloar also measured activities of the 
benzene-rubber system at different temperatures, 
thus enabling the separate evaluation of the 
partial molal entropy and partial molal heat of 
dilution. In dilute solutions the partial molal 
entropy is not in accord with Eq. (12), nor is 
the partial molal heat of dilution consistent 
with Eq. (16a) throughout the concentration 
range.*!°13 The former quantity approaches the 
theory as the concentration of rubber is in- 
creased. When these two somewhat erroneous 
equations are combined, however, a satisfactory 
free energy function is obtained, as Huggins has 
shown. 

The basis for this divergence between the 
present theory of polymer solutions and experi- 
mental entropies has become evident from 
further theoretical research into the nature of 
dilute polymer solutions.* It is unnecessary to 
go into the causes for the discrepancies further 
than to point out that in the theory as here (and 
previously) developed it is assumed, with refer- 
ence to the jth molecule to be added to the 
lattice, that all previously added molecules 
present a perfectly random distribution of cells 
occupied by polymer chain segments. This is not 
a satisfactory approximation in dilute solutions 
of high polymers. Here the occupied sites are 


4a The parameter yu in this paper corresponds to Huggins’ 
Mi. 
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present in loose ‘‘clusters,’’ each polymer mole- 
cule occupying a “‘cluster’’ of consecutively 
adjacent cells. Intervening regions between these 
“clusters” are completely vacant. 

For present purposes it is sufficient to point 
out that the theory, as given in previous papers 
and as extended here, yields an expression for the 
partial molal entropy which when combined with 
the » term gives a partial molal free energy of 
the solvent which agrees well with experiment. 
It is only necessary to modify the significance of 
uw to include, in addition to the heat term, con- 
tributions from other factors the origins of which 
are not yet clear.'* With yu redefined in this way, 
Eq. (17) appears to depict the partial molal free 
energy of the solvent satisfactorily. 

For the two-component system consisting of a 
single polymer species of size x dissolved in a 
solvent, the condition 


(AF ,/dn) = (dAF,/dN;) 


e  (@AP/a0:) = (x01/v2)(@AP 1/80.) 


necessarily applies. This connection permits 
derivation of the concentration dependent terms 
of AF, from AF. Thus, since Eq. (17) accounts 
satisfactorily for experimentally observed partial 
molal free energies of the solvent in two com- 
ponent systems composed of polymer and solvent 
(in which case v,=v2 and #,=<x), then it neces- 
sarily follows that 


AF,=RT[In v,+0,x(1—1/x) 
+px(1—v,z)?+const.] (18’) 


must be equally acceptable for the solute. 
Furthermore, Eq. (17) applies with approxi- 
mately the same value of u to all polymers of a 
given series regardless of molecular weight. 
Hence, Eq. (18’) must also be applicable to 
polymers of different molecular weights. 

These, then, are the experimental facts jus- 
tifying the use of Eq. (17) and, for solutions of 
one polymer component, Eq. (18’) as well. 
Returning to the theory, it will be recalled that 


4b Huggins’ treatment, which differs slightly from that 
of the writer, leads to a theoretically derived » which 
includes other terms in addition to the heat term. How- 
ever, in the opinion of the writer these are trivial in com- 
parison with the difference between the value of u required 
to_reproduce the observed activities and the value of 
AAl,/RT obtained from temperature coefficients. 


mg 
=e 
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no significant new assumptions are required in 
the generalization of the treatment to apply to 
heterogeneous polymers. It is justifiable, there- 
fore, to place equal confidence in the expression 


AF, =RT[In v2+02x(1—1/2,) 


+yx(1—v2)?+const.] (18”) 


for the partial molal free energy of x-mer' in a 
solution of a heterogeneous polymer. 

The terms of Eq. (18) which are replaced by 
an undetermined constant in Eq. (18’’) originate 
from the entropy of disorientation. The limita- 
tions of the theory in accounting for entropies in 
dilute solution have nothing to do with the 
validity of these terms in the free energy func- 
tion. This portion of the partial molal free energy 
expression (18) rests squarely upon the theory in 
its present state of development, and is neither 
supported nor contradicted by experimental 
results. 

The above arguments have been presented 
rather in detail in an effort to show clearly the 
justification for usings Eqs. (17) and (18), or 
(19), to represent the partial molal free energies, 
even in dilute solutions where the present theory 
fails to give correct entropies. In the following 
sections these equations will be applied.as semi- 
empirical approximations to the actual free 
energy functions. 


APPLICATION TO THE FRACTIONATION OF 
HIGH POLYMERS 


Frequently it is desirable to separate a high 
polymer into a series of fractions differing in 
average molecular weight. This is usually accom- 
plished by fractional precipitation. Either of two 
methods may be employed: Fractions of suc- 
cessively decreasing molecular weights may be 
obtained by adding a precipitant in suitably con- 
trolled increments to a dilute solution of the 
polymer maintained at constant temperature, 
the precipitate formed after each addition of 
precipitant being separated from the supernatant 
solution ; or a simple solvent medium is employed 
from which precipitation is induced by lowering 
the temperature in suitable stages. It is neces- 
sary, of course, to assure equilibration after each 
addition of precipitant, or after each decrease in 
temperature. In a few instances fractional solu- 


tion rather than fractional precipitation has been 
employed. 

In all of these methods, separation of suc- 
cessive fractions depends on a controlled altera- 
tion of uw in the critical region for precipitation. 
This critical region occurs when y» exceeds one- 
half.6 Except in the case of polymers which 
precipitate in crystalline form, the precipitates 
so obtained are liquid phases containing large 
amounts of solvent, although the concentration 
of polymer therein exceeds considerably the con- 
centration in the supernatant phase. Frac- 
tionation depends upon the partition of each 
polymer species between these two phases. 
Higher polymers preferentially occur in the more 
concentrated phase. It has been recognized that 
the resulting fractions are by no means sharp; 
each fraction contains a rather wide distribution 
of species, although it is less heterogeneous than 
the original polymer. The equations given above 
provide a means for calculating the degree of 
separation obtained in the fractionation of high 
polymers under various conditions. They may 
also be of assistance in clarifying the factors 
affecting the sharpness of separation. 

As indicated above, it is assumed in the fol- 
lowing analysis that no crystallinity, or related 
regular structure, occurs in the precipitated 
phase. Furthermore, all polymer species are 
assumed to possess the same composition and 
unit structure, their only differences residing in 
the lengths of their chains. This latter assump- 
tion assures that the same u value applies to all 
species. It will be assumed further that yp pos- 
sesses the same value in either phase. This last 
assumption in general will not be strictly true 
when the solvent medium comprises a mixture 
of solvent and precipitant. The two phases may 
contain different proportions of the precipitant. 
The magnitude of the differences in » for the two 
phases which may arise from this source has not 
been investigated. However, introduction of dif- 
ferent » values for the two phases would not 
change the general nature of the results derived 
below. When a single liquid is used for the 
solvent medium and the fractions are obtained 
by successive decreases in temperature, this 
assumption becomes unnecessary. 

For equilibrium between the two phases to 
exist it is necessary that the partial molal free 





energy of each constituent be the same in both 
phases. Using primes to signify quantities re- 
ferring to the precipitated phase, the conditions 
for equilibrium can be expressed as follows: 


AF,=AF\, AF,=AF/’. 
Substituting from Eqs. (17) and (19), or (18), 
In (1 —v2) +02(1 —1/%,) +uv2* 

=In (1—v2') +02'(1—1/2n')+uv2", (20) 
(In vz) /x+v2(1—1/Fn) +u(1—v2)? 
= (In 9’) /x-+09'(1—1/En')+u(1—v2')*. (21) 
Equation (21) can be written 


In (v,’/vz) = ax, (22) 
where 


a=v2(1—1/%,) —ve/(1—1/2,’) 
+ul(1—ve)?—(1—00')?]. (23) 


The quantity a can be regarded as a partition 
factor. The ratio of concentrations of a given 
species in the two phases varies exponentially 
with x. 

An equation equivalent to Eq. (22) has been 
obtained by Schulz.!® In deriving his relationship 
Schulz assumed that partition between the two 
phases is governed solely by differences in heats 
of solution of polymer molecules in the two 
phases. He assumed that this difference should 
increase in proportion to the molecular weight; 
thus the larger molecules should be preferentially 
distributed in favor of the precipitated phase 
where their heat of mixing with the medium will 
necessarily be lower. He neglected entropy effects 
such as have been taken into account by the 
present theory. Nevertheless, Schulz’s partition 
equation is equivalent to Eq. (22) above, and 
he has carried out significant calculations on the 
nature of the molecular weight distributions in 
fractionated polymers. The significance assigned 
by Schulz to the partition factor a must be 
revised, however. His definition of a@ in effect 
. includes only the u term of Eq. (23). 

The application of Eq. (23) to fractionation 
problems can be indicated as follows. Let V and 
V’ represent the volumes of the respective phases, 


16 G. V. Schulz, Zeits. f. physik. Chemie B46, 137 (1940); 
ibid. B47, 155 (1940). 
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and let Q, and Q,’ represent the absolute volumes 
of x-mer in them; thus Q,= Vv, and Q,’=V'v,’. 
Also, let Q,° equal the total volume of x-mer, i.e., 
Q0.°=Q.+Q,’, and let the ratio of volumes of 
precipitated and supernatant phases be repre- 
sented by r= V’/V. Then according to Eq. (22), 


Q.=Q.°/(1+re**), (24) 
Q.,' =re**Q,°/(1+7re**). (24’) 


Whereas the relative concentration of each 
species in the two phases is governed by a, their 
relative amounts depend also on r. 

The conditions of fractionation depend jointly 
on a@ and r. Both the location of the mean 
molecular weight about which the polymer is 
divided between the two phases and the sharp- 
ness of the separation depend on the particular 
combination of these parameters. Let us charac- 
terize the location of the region of separation by 
the molecular weight, or x value, at which 
Q.=Q,'. Let this particular value of x be desig- 
nated by xa. Then according to Eqs. (24) and 
(24’), 

exp (ax,) =1/r. 


Further, as an index of the efficiency of the 
separation let us consider the slope of the curve 
representing the proportion of x-mer in the 
precipitated phase plotted against x, at the 
point x=x,. According to Eq. (24’) this quantity 
is given by 


[d(Q.'/Q:°)/dx ],,= 0/4. 


The relative efficiency ¢« can be defined as the 
increase in Q,’/Q,° per unit relative change in x 
when x=X, 


e=([d(Q.'/Q:°)/d In #1. 
€=ax,/4= (1/4) In (1/r). (25) 


Thus the relative efficiency of separation in- 
creases as the log of the reciprocal of 7, and is 
independent of the molecular weight region in 
which separation occurs.'5 

From Eq. (25) it is clear that the volume of 
the supernatant phase relative to the volume of 
the precipitated phase should be made as great 
as possible. The importance of this point gains 
further emphasis from the fact that, according 
to Eqs. (24) and (24’), the concentration of 
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every species including the lowest is greater in 
the precipitated phase than in the more dilute 
phase, although for components well below the 
molecular weight region of separation these con- 
centrations will be nearly equal. Efficient elimina- 
tion of lower polymers from the precipitated 
phase is possible only in the presence of a very 
large supernatant phase. Consequently, the solu- 
tion from which precipitation takes place should 
be very dilute. 

From a consideration of the phase equilibrium 
diagram (Fig. 1) depicting the concentrations of 
the phases in equilibrium as a function of » which 
in turn depends on temperature or amount of 
precipitant, it is obvious that the more dilute the 
solution the more concentrated will be the pre- 
cipitated phase. Thus, dilution not only increases 
the volume of the supernatant phase, but it also 
produces a ‘“‘deswelling’” of the precipitated 
phase. Both of these changes contribute toward 
the desired decrease in r. 

The phase equilibrium diagrams for high- 
polymer solutions are extremely unsymmetrical 
(Fig. 1). The critical concentration at which 
partial miscibility first develops with increasing 
uw occurs at a low concentration of polymer. This 
critical concentration for a homogeneous polymer 
is given® by 


Vaceritical) = 1/(1++/x). 


It moves in the direction of more dilute solutions 
as the molecular weight (x) is increased. From 
an examination of the shift with molecular weight 
in the phase equilibrium curves shown in Fig. 1 
for homogeneous polymers, it is clear that the 
greater the molecular weight, the greater the 
degree of swelling of the precipitated phase in 
equilibrium with a supernatant phase of given 
concentration. Therefore, with increase in the 
molecular weight at which precipitation occurs, 
it becomes increasingly important to use dilute 
solutions. In fact, if the concentration is not much 
less than ve (critical) as given above, very little 
fractionation is to be expected. 

The above discussion indicates the nature of 
the dependence of fractionation characteristics 
on the parameters a and 7. Unfortunately, these 
are not simple functions of externally controlled 
variables such as », which depends on the amount 
of precipitant or temperature. Equation (23) 


gives an indication of the complexity of the de- 
pendence of a on the concentrations of polymer 
in the two phases, on the number average molecu- 
lar weights and on u; the quantities other than yz 
on the right of Eq. (23) depend on the initial dis- 
tribution and on the dilution. The following dis- 
cussion is aimed at a simplification of these 
inter-relationships. 

The ratio of the volumes of polymer (all species 
included) in the two phases, >> Q.’/>> Q:, must 

x x 


equal rv’/v2, from which it follows that 
DX [Q2°e**/(1+re*”)] v2" 


=—=y, 26 
5 (0./(1+re*)] ais 


where the symbol y has been introduced to repre- 
sent this ratio. Equation (26) relates r, a, and the 
initial polymer distribution (the V,°’s) on the 
one hand with the relative total polymer concen- 
tration v2'/v2 on the other. 

The number average molecular weights are 
similarly related to these three quantities. 


En= Lo v2/L (vz/x) =02/L0 (vz/x). 


The quantity v2/Z, which appears in the equation 
for a, and in subsequent equations is given by 


v2/En=), (vz/x) 
=(1/V) © (Q.°/x(1+re%)]. (27) 





Similarly 
v2! /En'=(1/V") D [(Qzre%*)/x(1+ree7) J. (28) 


Replacing V’ with rV and subtracting Eq. (27) 
from (28), 


(v9 /E,’) — (v2/En) 
=(1/V) ¥ [Q.%(e**—1)/x(1+re2*) ]. 


It will be convenient to define a quantity 6 as 
follows: 


5=(1/v2)[(v2" /En’) — (v2/Fn) J. 


Since 
V=D (Qz/v2) =(1/v2) DY [Q2°/(1+re2*) ], 


it follows that 6 is given by the following function 
of the initial distribution, 7 and a: 


© [0.%(er*—1)/x(1+r0**)] 
x (0.°/(i+re2)] 





(29) 
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Returning now to Eq. (23) and solving for y, 


= a+ (v2" — V2) — 62 (30) 
2 (v2! —v2)[2— (v2 +02)] 


Substituting this expression for » into Eq. (20), 
which as yet has not been made use of, and 
rearranging, 


[1 —(v2’ +02) /2] In [(1 —v2)/(1 —v2")] 
= (9' —v2) + a(v2' +02) /2 — v2. 





Series expansion of the logarithmic term fol- 
lowed by combination of terms resulting from 
multiplication by the factor in brackets ulti- 
mately leads to 


Cvs" | 1 ae 
V2 —v —— +} 
i, <i e 

3v9/2? +4 v0/v2+ 3022 

4-5 

(n+1)v9!"+ 2nv2'"—v2+ - - + (n+1)v2”" 

+ +] 
(n+2)(n+3) 
= ave’ +02) /2 — dv. (31) 








The infinite series in braces is very nearly equal 
to 1/6(1—v2’—v2) so long as v2 and v2’ are small. 
Making this substitution in Eq. (31), 


(v0" —v)?/12(1 —v~! —V2) = a(ve’ +02) /2 - 6V2. (32) 


The percentage error introduced by the ap- 
proximation will be of the order of (v2’?/10)100, 
since v2<v.’. Letting the ratio of the polymer 
concentrations in the two phases be represented 
by y=v2'/v2, Eq. (32) may be expressed alter- 
natively as 


(y—1) 902? = 12[1 —v2(y+1) J 
XO(v+1)e/2—6]. (32’) 


The external variables associated with separa- 
tion of heterogeneous polymers by fractional 
precipitation are the initial polymer distribution, 
the value of u (determined by temperature and 
amount of precipitant), and the average polymer 
concentration, or concentration of polymer prior 
to separation into two phases. This average con- 
centration is given by (v2+7v2’)/(1+r). The com- 
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plexity of the relationships derived above is such 
that it has not been possible to deduce a and r 
directly from these independent variables. How- 
ever, it is possible to obtain solutions by assigning 
arbitrary values to the dependent variables r 
and a and then solving for the physically inde- 
pendent variables. Given the initial distribution 
(the V,°’s) and assumed values of r and a, 
v2’ /v2=~ can be obtained by executing the sum- 
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Fic. 1. Equilibrium phase compositions for solvent: 
homogeneous polymer mixtures in the region of partial 
miscibility for several values of x (relative molecular 
weight). Curves calculated from Eqs. (33) and (30). The 
pair of solutions for a given value of u represent the cal- 
culated compositions of the phases in equilibrium. Since u 
increases with decrease in temperature, the ordinate may 
be regarded as a temperature scale, and the curves should 
resemble temperature-composition diagrams for the binary 
liquid system. 


mations of Eq. (26). Similarly, 6 may be obtained 
from Eq. (29). Substituting these values of y 
and 6 into Eq. (32’), the resulting quadratic 
equation can be solved for ve. The value of v2’ 
is given by yve. The average concentration can 
then be deduced, and y» can be computed from 
Eq. (30). This completes the deduction of the 
independent variables from the assumed values 
of r and a. Since the polymer distribution in each 
phase is governed, according to Eqs. (24) and 
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(25), by r and @ and the initial distribution, the 
solution for this particular set of independent 
variables, the V,°’s, u, and the average concen- 
tration, is complete. 

Equation (32), or (32), is useful also in deriv- 
ing phase equilibria for homogeneous polymer- 
solvent systems. When only one polymer 
species composed of x segments is present, 
ve=v, and v2’ =v,’. Hence 


a=(In y)/x, 6=(y—1)/x 


and Eq. (32’) reduces to 
(y—1)*02? = 12[1 —v2(y +1) ] 
XL(y+1)(In y)/2—(y—-1) ]/x. (33) 


For any given value of x, this equation may be 
employed to obtain the concentrations of the 
two phases which are in equilibrium with one 
another. By substituting an arbitrary value for 
the ratio y of these concentrations, solution of 
the quadratic equation (33) gives the concentra- 
tion v2 of the more dilute phase. Multiplication of 
ve by y gives, of course, the concentration in the 
more concentrated phase. The corresponding 
value of » can be found from Eq. (30). This pro- 
cedure, based on the approximate equation (33), 
which is a’ modification of (32), is very much 
simpler than the method of successive approxi- 
mations previously employed.® 

Solutions of Eq. (33) are shown in Fig. 1 
plotted against yu, for several values of x. The or- 
dinate scale is plotted in the direction of de- 
creasing yw. Since uw decreases with increase in 
temperature, these plots should resemble tem- 
perature-composition diagrams for the binary 
liquid systems. The positive direction of the 
ordinate axis corresponds to an increase in tem- 
perature, or to a decrease in the proportion of 
precipitant at constant temperature. 


CHEMICAL EQUILIBRIA IN POLYMERIC 
SYSTEMS 


An appropriate formulation of an equilibrium 
constant for a chemical reaction involving poly- 
mer formation, condensation, or degradation can 
be secured with certainty only when the de- 
pendences of the partial molal free energies (or 
the activities) of the reacting components on 
their concentrations are known. This is obvious 


from reference to the thermodynamic origins of 
the equilibrium constant. In this section the 
application of the above thermodynamic func- 
tions to the equilibria governing associative and 
degradative reactions of polymer molecules will 
be discussed. 

Various condensation polymers which may be 
formed and degraded reversibly serve to illus- 
trate the type of equilibria considered here. 
Polyesters, for example, are synthesized from 
their monomers by the reversible esterification 
reaction 


xHO—R—COOH@HO[—R—COO— ].H 
+(x—1)H,0. 


The polymer molecules actually are not formed 
kinetically by the simultaneous condensation of 
x monomers; they are formed in an irregular 
stepwise succession of condensations, each of 
which is reversible,!® 7 


HO[—RCOO],—H +HO[—RCOO—],H2 
HO[—RCOO}],4,—H+H.0. 


Reactions of this character can be represented 
by the scheme'® 


A—B@A—BA—BA—etc., 


where A and B are the functional groups which 
co-react to form the inter-unit linkage AB. In 
this process a molecule such as water may or 
may not be split out. The following treatment is 
also applicable, with minor modifications to the 
polymerization of a pair of bifunctional monomers 
A—A and B—B, where the complementary 
functional groups reside in different monomers. 
Vinyl type polymerizations could be included, 
although these processes generally are not 
reversible in the thermodynamic sense. 

As the standard reference state for each 
molecular species, a perfectly ordered hypo- 
thetical crystal will be chosen in which the 
molecules are arranged in perfect order. The 
standard states for the various species will 
possess similar structures insofar as the arrange- 
ment and packing of the units are concerned. 


1% P. J. Flory, J. Am. Chem. Soc. 61, 3334 (1939). 
17P, J. Flory, J. Am. Chem. Soc. 64, 2205 (1942). 
18 P. J. Flory, J. Am. Chem. Soc. 58, 1877 (1936). 











Furthermore, the packing density will be so 
chosen that the heat change in passing to the 
actual liquid state with randomly arranged 
molecules is zero. It should be remarked that 
these conditions are introduced as definitions of 
an arbitrary standard state, and not as assump- 
tions. 

We now assume that with the standard states 
so defined the standard free energy change in a 
reaction 


x-mer-+y-mer=(x+y)-mer (34) 


is always the same regardless of x and y, and 
equal to AF®. In other words, it is assumed that 
the standard free energy change for a reaction 
leading to the formation of one AB linkage (e.g., 
one ester group) is invariably equal to AF°. This 
assumption can be expected to apply at all 
values of x and y when the functional groups A 
and B are separated sufficiently in the monomers 
so that the reactivity of A is influenced neg- 
ligibly by previous reaction of the B group of 
the same unit, and vice versa. Furthermore, the 
intermolecular attractions of the A and B groups 
should not differ too much from that of the AB 
linking group. Even when these conditions are 
not met, the above assumption can be expected 
to apply satisfactorily except for low values of 
x and y. 

The free energy of an x-mer admixed with 
other species, and possibly with a solvent as 
well, is given by 


F,= F,°+AF,, 


where F,’ is its free energy in the standard state 
and AF, is the free energy change in transporting 
one mole of x-mer from the standard state to an 
infinite amount of the mixture, or solution. In 
other words, AF, is the partial molal free energy 
of the x-mer in the mixture with respect to the 
standard state, as given by Eq. (18). According 
to the definition of the standard states and the 
assumptions introduced above, 


FP=xFP+(x—1)AP, (35) 


where F;,° is the free energy of the monomer in 
the standard state and AF°, as previously defined, 
is the standard free energy change per mole of 
inter-unit linkages AB in the molecule. The 
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standard states have been defined with the 
purpose of view of obtaining an equation of the 
form of Eq. (35) wherein AF°® is a constant. 
Hence, 


F,=x(FY+AF) —AF°+A4P,. 


The total free energy of the mixture exclusive 
of inert solvent is given by 


F=>, N.F,=(FY+AF*) > xN,z 
When the polymerization (or degradation) reac- 
tions have reached a state of equilibrium, dF=0. 


Taking the differential of Eq. (36) and equating 


to zero, 7 
> (AF,—AF*)dN,=0, (37) 


since constancy of the total amount of material 


requires that 


and since > N,d(AF,)=0, according to the 
Gibbs-Duhem equation. Equation (37) ex- 
presses the necessary condition for equilibrium. 

If the process concerned is a simple one such 
as is expressed by Eq. (34), where only three 
species are involved, then 


dN,=dN,= —dNr4y, 


and for the equilibrium condition we have, ac- 
cording to Eq. (37), 


AF,,,—AF,—AF,+AF°=0. 


After substituting from Eq. (18) for these partial 
molal free energies there is obtained 


In K=—AF°/RT+1n [(y—-1)o], (39) 


Urt+y 

<A +#* 

Vay/ \x+y 
Alternatively 


In K’= —AF°/RT+In [(y—1)0 Vi], (39) 


where 





where V; is the molar volume per structural unit 
A—B, and 
K! =Crpy/Caly, (40’) 


the c’s representing the respective concentra- 
tions in moles per unit volume. Inasmuch as y 
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is a function of the liquid structure only and o 
and V, are constants characteristic of the types 
of molecules concerned, evidently either K or K’ 
should be an acceptable equilibrium constant. 

It should be noted that use of mole fractions 
instead of molar concentrations in expressing the 
equilibrium constant (except, of course, when the 
equilibrium constant is dimensionless) is irrecon- 
cilable with the statistical mechanics of these 
mixtures.’® In the absence of diluent, conversion 
to mole fractions X, through the substitutions 


v,=X;z >, N,z/> xNz, 
etc., in Eq. (40) yields 
[X24/X2Xy]=K Dd N:/L x«N.=K > Nz/No, 


where the total number of polymer segments, 
> xNz, is replaced by No. Since >> Nz varies 
with the degree of displacement of the reaction, 
the quantity in brackets should not be employed 
as an equilibrium constant. 

Proceeding now to the more complex equi- 
librium between all possible species M,, or 
(A—B),, where all values of x are considered, 
we will again employ the equilibrium condition 
as expressed by Eq. (37). According to Eqs. (18) 
and (38), 


> APAN,=RT ¥ {In (Nzv2/No) 
+In [(y—1)o]}dN.. (41) 
Substituting Eq. (41) in (37) 


Xj In (N2v2/No)+in [(y—1)o] 


AF°® 
—— ;dN,=0. (42) 
RT 
Solving Eqs. (38) and (42) by the Lagrange 
multiplier method, 
N.z=(No/v2)e%e?, (43) 
where X is the multiplier and 


B=In [(y—1)a]—AF*/RT. 


19R. E. Powell and H. Eyring, J. Am. Chem. Soc. 65, 
648 (1943), in their analysis of the equilibria in liquid 
sulfur between ring (Ss) and chain forms have expressed 
equilibrium constants in terms of mole fractions. According 
to the present treatment their equilibrium “constants” 
are not to be regarded as constants. 


It is found upon substituting Eq. (43) into 
the equation }> xN,= No, and carrying out the 
summation over integral values of x from 2 to 
o, that 

e® =e /vo(1—e)?. (44) 


Also, the total number of molecules >> N, 
must equal N,o(i—/p), where p is the “extent of 
reaction,”’ or fraction of the initial A, or B, 
groups converted to AB. According to Eq. (43), 
after substituting Eq. (44) for e’, 


o Ne=Ni(i—e>). 


Hence, e~* may be identified with the extent of 
reaction p. Making this substitution and intro- 
ducing Eq. (44) in (43), 


N.z=No(1—p)*p7", (45) 


which is the equation, previously derived’: !* by 
other methods, for the molecular size distribution 
in linear condensation polymers. 

Tobolsky*® has derived the above molecular 
size distribution equation by a similar procedure 
based on the same thermodynamic equations.® 
He sought to introduce the effect of heat of 
polymerization into the equations, with particu- 
lar reference to vinyl type polymerizations which 
invariably are exothermic. In so doing he errone- 
ously employed the heat of reaction where a 
standard state free energy change obviously 
should be used. Consequently, an equilibrium 
constant derived from his treatment would de- 
pend on the former, rather than on the latter. 
The error can be corrected by substituting AF° 
for P (heat of reaction) in Tobolsky’s equations. 
If pursued further, his treatment corrected in 
this manner would lead to conclusions equivalent 
to those reached here. 

According to Eq. (44), 


B=In [o(y—1) ]—AF°/RT=In K, (46) 
where 


K= 2 = 2. 
ee p/v2(1—p) (47) 


In K’=In [o(y—-1)Vi]—AF*/RT, (46’) 


where 
K’'=pV/v2(1—p)?=Cap/Cace. (47’) 


Cap, Ca, and Cz represent the concentrations of the 
respective functional groups in moles per unit 


20 A. V. Tobolsky, J. Chem. Phys. 12, 402 (1944). 
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volume. Thus, polyfunctional condensations may 
be treated as reactions between functional groups 
without regard for the particular molecular 
species involved, a point of view set forth in 
previous publications.!*'8 The entropies of mix- 
ing of long chain polymers are such as to require 
no revision in this postulate as applied to chem- 
ical equilibria between polymeric reactants.”° 

It.can be shown further that even the mag- 
nitude of the equilibrium constant should not be 
altered by the polymeric nature per se of the 
components in such equilibria. That is to say, 
the percentage esterification of a glycoldibasic 
acid mixture at equilibrium should be about the 
same as is found between analogous monoalcohols 
and carboxylic acids under equivalent conditions 
of temperature and dilution. A monofunctional 
reaction of this type comes within the scope of 
the treatment of the reaction (34). If there are 
initially co moles of each of the two reactants per 
unit volume, then 


Ca=Ce=C(1—p), Cap=Cop 
Co=V2/( Vat Vz), 


where V4 and Vz are the molar volumes of the 
reactants. Substituting in Eq. (40’) 


K'=p/(1—p)*co= (Vat Va) /v2(1—p)?, 


which is identical with Eq. (47’) for the bifunc- 
tional reaction, except for the obvious replace- 
ment of Vi by (Va+ Vs). 

In the above discussion the role of a by-product 
of the polymerization reaction, such as the water 
formed during polyesterification, has been dis- 
regarded. It is obvious from the above discussion, 
however, that equilibrium constants for con- 


and 
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densation processes from which a by-product is 
released may be written in an analogous manner, 
the concentration of the by-product appearing in 
the numerator. 

These generalizations concerning equilibria in 
polymer reactions can be applied to the fre- 
quently encountered reversible ring-chain equi- 
libria. For example, the six-membered lactide 
ring” can be reversibly converted to a linear 
polyester 


co 
—, 
CH;—CH ’ 
‘ 
O CH—CH; 


CO 


iu 


[Bod 
HO—L—CH—CO—OJ,—H. 


In this system two equilibria are to be con- 
sidered: the ring-chain equilibria indicated 
above, and the polyesterification equilibrium, 
involving water as a by-product, which deter- 
mines the average length of the chains. These 
chains will possess the size distribution given 
above, and this phase of the equilibrium may be 
treated as a reaction between functional groups. 
Other examples include numerous six-membered 
lactone rings which can be reversibly converted 
to chain polymers, the (PNCl:)3—(PNCl2)- 
reaction,” aldehyde polymers, cyclic lactams, and 
the ring-chain sulfur transformation already 
referred to.’® 


21'W. H. Carothers, G. L. Dorough, and F. J. van Natta, 
J. Am. Chem. Soc. 54, 761 (1932). 
2 See K. H. Meyer, reference 4, p. 52. 
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The Dielectric Properties of g-Lactoglobulin in Aqueous Glycine Solutions and in the 
Liquid Crystalline State 


T. M. Suaw, E. F. JANSEN, AND H. LINEWEAVER 
Western Regional Research Laboratory,* Albany, California 


) (Received July 17, 1944) 


- The dielectric constant e’ and loss factor e” of solutions of §-lactoglobulin dissolved in 
aqueous 0.48, 1.5, and 2.5 molar glycine, and the dielectric constant of 8-lactoglobulin liquid 
' crystals have been measured for frequencies between 0.01 and 5.0 megacycles. The data show 


a broader frequency range of dispersion and absorption and a smaller value of ies than is 
predicted by the Debye theory for a system characterized by a single relaxation time. The 
disparity between the experimental data and the Debye theory increased with the concentra- 
tion of lactoglobulin and was greater than expected on the basis of the small effects of protein 
concentration on other physical properties of B-lactoglobulin, e.g., sedimentation constant. 
Dielectric constants for “‘zero”’ and “‘infinite” frequencies have been obtained by the Cole cir- 
cular-arc method of extrapolation. The frequency dependence of e’ and e’’ was found to be well 
represented by the empirical equation of Cole and Cole, e’ —ie’’ =€..+ (€o—€x)/[1+(iwro)**]. 
In this equation €9 and ¢, are, respectively, the dielectric constants for zero and infinite fre- 
quencies, w= 27 times the frequency, 70 is a generalized relaxation time, and a is a constant, 
0<a<1. For a=0, the equation reduces to that of Debye for a system characterized by a 





l. single relaxation time. For the solutions of 8-lactoglobulin a was found to range from 0.08 to 
0.4, depending upon the lactoglobulin and glycine concentrations. The relaxation time varied 

1- from 7.5 to 25X10- sec. for a change in lactoglobulin concentration of 3.5 to 94.0 grams per 

d liter. The dipole moment is estimated to be 790+26X 10~ e.s.u. 

1, 

[- 

3€ INTRODUCTION 


of 0.48, 1.5, and 2.5 molar glycine containing up 


n HE importance of electrical properties for to about 10 percent of lactoglobulin and for 
ye the characterization of proteins and other lactoglobulin liquid crystals containing about 33 
" macromolecules has been established by many  Pércent protein. 
- experimental studies. Exploratory experiments ew 
directed toward the characterization, by dielec- 
. tric constant measurements, of such molecules Materials Studied and Their Preparation 
l . . . 
and their enzymic degradation products, have . 7 al — 
ly led eff y oe ag : The £-lactoglobulin used in this investigation 
revealed eifects at high protein concentrations ; , 
not previously reported. In view of the im- was prepared from cow's milk by ie maies 
, . described by Palmer.? The preparation was 
portance of these effects for the interpretation of ” ; ee ; 
- dielectric constants in terms of molecular size crystallized five times by dialysis. Stock solutions 


shape, and charge distribution, they have been of lactoglobulin in glycine were prepared by 


systematically investigated for the protein B- ‘ 
lactoglobulin. Extensive dielectric measurements © ; W= I/T. 
of dilute solutions containing one percent or less ; 

of this protein have been reported by Ferry and 
Oncley.! In the present paper dielectric constant 
measurements are reported for aqueous solutions 








* One of the laboratories of the Bureau of Agricultural a 
and Industrial Chemistry, Agricultural Research Adminis- 


; 
bs 
' 
k 


tration, U. S. Department of Agriculture. The work re- 
ported was begun in 1939 in the Enzyme Laboratory of 
the Bureau and was completed in the Western Regional 
Research Laboratory. 

1J. D. Ferry and J. L. Oncley, J. Am. Chem. Soc. 63, 
272 (1941). 


Fic. 1. Cole circle diagram for complex dielectric con- 
stant e* =¢’—ie’’. The case illustrated is for a>0. For a=0 
the center of the arc falls on the e’ axis. The maximum 
value of e”’ occurs for w=1/ro. 


2A. H. Palmer, J. Biol. Chem. 104, 359 (1934). 
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Fic. 2. Dielectric constant e’, and loss factor e’’, for solutions of B-lactoglobulin 
in 0.48 molar aqueous glycine at 20.0°C. Data are shown for lactoglobulin con- 
centrations C equal to 3.48, 10.4, and 20.9 grams per liter. 


adding crystals of lactoglobulin to glycine solu- 
tions. After analytical determination of the pro- 
tein and glycine concentrations, the glycine con- 
centration was adjusted to the desired value by 
the addition of solid glycine. To obtain solutions 
containing different amounts of lactoglobulin the 
stock solutions were diluted with glycine solu- 
tions containing the same amount of glycine as 
the stock solutions. In this way solutions con- 
taining various amounts of lactoglobulin dis- 
solved in 0.48, 1.5, and 2.5 molar glycine were 
obtained. The concentrations of glycine and 
lactoglobulin were calculated from determina- 
tions of total nitrogen and protein nitrogen in the 
solutions. The pH of the solutions was between 
5.6 and 6.0. 

The liquid crystals* were prepared by dialyzing 
concentrated solutions of lactoglobulin containing 
sodium chloride until sufficient salt was removed 
to cause the protein to separate as a new phase. 
Excess water was removed from the protein phase 
by centrifugation and the liquid crystals were 
transferred to the dielectric cell. During the 
entire process of preparation, the temperature 
was maintained between 10° and 15°C. 


Electrical Measurements 


A detailed description of the apparatus and 
methods used to determine dielectric constants 
has been published previously.’ 


* This material is the ‘‘oil’’ described by Palmer (refer- 


ence 2). 
3T. M. Shaw, J. Chem. Phys. 10, 609 (1942). 


Briefly, dielectric constants and dielectric loss 
factors were obtained from substitution measure- 
ments of the capacitance and conductance of a 
micrometer-type cell containing the solution of 
lactoglobulin. Measurements were made for about 
twenty frequencies spaced more or less uniformly 
over the frequency range of 0.01 to 5.0 megacycles. 
The bridge equipment employed for these meas- 
urements is similar to that used by Oncley.‘ The 
test cell is a modified form of the electrolytic cells 
used by Fricke and Curtis’ and by Cole and 
Curtis.° 

In the capacitance measurements, lead in- 
ductance and electrode polarization were the 
principal sources of error and appropriate cor- 
rection measures, already described in detail,* 
were taken. From calibration measurements in 
which glycine solutions of known dielectric con- 
stant were employed, the probable error in the 
dielectric constants was estimated to be +0.5 
unit for the protein solutions and about +1.0 
unit for the liquid crystals. 

In the conductance measurements an empirical 
correction was made for the error introduced at 
the higher frequencies by lead inductance and 
skin effect in the wire-wound resistance stand- 
ards. These corrections were obtained from a 
series of calibration measurements in which the 
dielectric test cell* filled with aqueous KCI solu- 


4J. L. Oncley, J. Am. Chem. Soc. 60, 1115 (1938). 
( & 5 Fricke and H. J. Curtis, J. Phys. Chem. 41, 729 
1937). 

®K. S. Cole and H. J. Curtis, Rev. Sci. Inst. 8, 333 
(1937). 
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Fic. 3. Dielectric constant ¢’, and loss factor e’’, for solutions of 6-lactoglobulin in 1.5 molar aqueous glycine at 20.0°C. 
Data are shown for lactoglobulin concentrations C equal to 12.7, 21.2, 38.2, and 63.7 grams per liter. 


tions was used as a standard electrolytic resistor. 
The conductance of the electrolytic resistor was 
determined accurately at low frequencies by 
bridge measurements. It was assumed that the 
change in conductance of the electrolytic resistor 
over the frequency range of 0.01 to 5.0 Mc is 
negligible. Experimental tests’ of the predictions 
of the Debye and Falkenhagen theory of elec- 
trolytic conductance indicate this assumption to 
be valid within the precision of the conductance 
measurements here reported. The over-all prob- 
able error in the conductance determinations was 
estimated to be +0.5 percent. 

For the purposes of the present work, the 
dielectric loss factor e’’ is defined by the equation, 


é’’ = (G—Gpo)d/Ka, (1) 


where G is the conductance at angular frequency 
w, Go is the limiting conductance at low fre- 
quencies, K is a constant for the dielectric test cell 
equal to 0.1700+0.0005 yuyuf cm*, and d is the 
spacing of the test cell electrodes. 

The use of Eq. (1) is based upon the assumption 
that the a.c. conductance of the dielectric ap- 
proaches a constant limiting value Go at low 
frequencies removed from the frequency region ‘of 
dielectric absorption. Values of Go were obtained 
in the present work from graphs of G versus log w. 


70. M. Arnold and J. W. Williams, J. Am. Chem. Soc. 
58, 2616 (1936). 





The rate of change of G with frequency between 
10‘ and 2 X10# cycles, in all cases except for the 
liquid crystals, did not exceed the probable error 
in the conductance measurements (0.5 percent). 
It should be noted from Eq. (1) that even small 
errors in Go may lead to very large errors in e’”, 
particularly at the lower frequencies or for solu- 
tions where the total energy absorption is small 
and G is nearly equal to Gp. 


RESULTS AND DISCUSSION 


A convenient quantitative description of the 
dispersion and absorption of many dielectrics is 
afforded by the empirical equation recently 
described by Cole and Cole,’ 


(e* —€x)/(€o—€x) =[1+(twro)'@ J}. (2) 


In this equation ¢* is the complex dielectric 
constant equal to e’ —ie’’, where ¢’ is the ordinary 
dielectric constant, ¢’’ is the a.c. loss factor, w is 
the angular frequency equal to 27 times the 
electrical frequency »v, €9 and ¢ are the limiting 
values of e’ for zero and infinite frequencies, re- 
spectively, 7o is a generalized relaxation time, and 
a is a parameter ranging from 0 to 1. For a=0, 
Eq. (2) reduces to the equation given by Debye’ 
for polar dielectrics. 

8K. S. Cole and R. H. Cole, J. Chem. Phys. 9, 341 
(1941). 


®P. Debye, Polar Molecules (Chemical Catalog Com- 
pany, New York, 1929). 
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Fic. 4. Dielectric constant e’, and loss factor e’’, for solutions of B-lactoglobulin in 2.5 molar aqueous glycine at 20.0°C. 
Data are shown for lactoglobulin concentrations C equal to 11.8, 23.5, 47.0, 70.5, and 94.0 grams per liter. 


The real and imaginary components of Eq. (2) 
are 


(€’ —€x)/(€0—€x) =${1—sinh (1—a)x/ 
[cosh (1—a)x+sin (ar/2)]}, (3) 
and 


e’’/(€9—€x) = 3 cos (am/2)/ 
[cosh (1—a)x+sin (am/2)], (4) 
where 
x=log. (w7o). 


The principal advantage of Eq. (2) is that it 
provides a quantitative description of the experi- 
mental behavior of dielectrics. As is well known, 
the classical Debye equations? fail to provide, in 
many cases, more than a qualitative description 
of the behavior of the dielectric constant and loss 
factor.® 

Experimental data which satisfy Eq. (2) yield 
a circular arc when e” is plotted against e’. If a 
semicircle is obtained with center on the ¢’ axis 
then a=0 and, as previously noted, Eq. (2) re- 
duces to the Debye equation for a polar dielectric 
with a single relaxation time. If, however, an arc 
with center depressed below the e¢’ axis is ob- 
tained, then a>0 and a distribution of relaxation 
times is indicated. This situation is illustrated in 
Fig. 1 where a graphical method of obtaining a is 
represented. 


Plots of e«’’ versus e’ for solutions of B-lacto- 
globulin in aqueous glycine are shown in Figs. 2-4. 
In all cases the data are well represented by 
circular arcs with centers depressed below the ¢’ 
axis. It is immediately apparent that the angle by 
which the center of the arc is depressed below the 
e’ axis varies with the composition of the solution. 
To investigate this dependence, values of a have 
been obtained for each of the circle diagrams in 
the manner outlined in Fig. 1. The resulting data 
are presented in Fig. 5 where a is plotted as a 
function of the concentrations of 8-lactoglobulin 
and glycine. Figure 5 shows that for all protein 
concentrations studied a>0 and hence the ex- 
perimental data cannot be represented by the 
Debye equations which apply only to a system 
with a single relaxation time. Moreover, a in- 
creases with increasing protein concentration and 
with decreasing concentration of glycine in the 
solvent. The extent of the failure of the Debye 
equations to represent the experimental data 
becomes more apparent if use is made of Eqs. (3) 
and (4). In Figs. 6 and 7, typical dispersion and 
loss factor curves are shown for solutions of 
8-lactoglobulin representing extremes of concen- 
tration. The solid line curves were calculated 
from Eqs. (3) and (4) using the parameters €0, €., 
and a obtained from the circle diagrams of 





)- 





Fig. 3.!° The broken line curves marked a=0 are 
the so-called Debye curves. The agreement be- 
tween experimental data and the curves calcu- 
lated from Eqs. (3) and (4) for a>0 is good in all 
cases. 

Dielectrics characterized by a distribution of 
relaxation times!” yield circle diagrams with 
a>0O, similar to those shown in Figs. 2-4 for 
8-lactoglobulin solutions. 

In all cases such systems can be characterized 
by a symmetrical distribution of relaxation times 
about a central or most probable relaxation time 
to. The numerical value of 7o can be obtained 
from the experimental data in a number of 
ways. Cole and Cole suggest a graphical method 
based on Eq. (2).* Probably the most con- 
venient method is that described by Fuoss and 
Kirkwood." The relaxation time is calculated 


from the frequency »», for which cosh! (€max/e’’)=0 
according to the relation t>=1/(27vm). In the 
present work the frequency v», was determined 
graphically. In the few instances where a com- 
parison of the Cole method with that of Fuoss 
and Kirkwood was made, no significant differ- 
ences in 7) were found. 

The relaxation times obtained for the B-lacto- 
globulin solutions are shown in Fig. 8 where 7 has 
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Fic. 5. Variation of « with mole fraction ratio R. 
(R equals mole fraction 8-lactoglobulin divided by mole 
fraction of glycine.) 


©Since the fourth parameter 7» needed in the calcu- 
lation could not be obtained precisely from the e’’ versus ¢’ 
plots, the required values were obtained from graphs of 


” 

cosh (€max/e’’) versus log v, according to a method of Fuoss 
and Kirkwood, described below. 

1R, M. Fuoss, J. Am. Chem. Soc. 63, 2401 (1941). 

121t should perhaps be noted that the e’ and e’’ data of 
W. R. Conner, J. Chem. Phys. 9, 591 (1941), for a poly- 
disperse solution of maple lignin in chloroform when 
plotted according to the Cole method yield a circular arc 
with a=0.25. 

* See reference 8, Eq. (6). 

13R. M. Fuoss and J. Kirkwood, J. Am. Chem. Soc. 63, 
385 (1941). 
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Fic. 6. Frequency variation of dielectric constant and 
loss factor, according to Eqs. (3) and (4), for solution 
containing 11.8 grams per liter of B-lactoglobulin dissolved 
in 2.5 molar aqueous glycine at 20.0°C. 


been plotted as a function of the 6-lactoglobulin 
concentration C. A correction for the solvent 
viscosity has been applied to 7» to reduce the 
relaxation times to a comparable basis. According 
to Debye the relaxation time varies with the 
viscosity and absolute temperature of the solvent. 
Accordingly the values of 7» have been corrected 
to correspond to the viscosity of water at 20°C by 
the equation, 
T=T0(Nw/Ns), 


where 7 is the relaxation time in water at 20°C, 
Nw = viscosity of water at 20°C, »,=viscosity of 
solvent corresponding to the observed relaxation 
time 79. 

The data in Fig. 8 show that 7 increases with 
increasing concentration of B-lactoglobulin but is 
independent of the glycine concentration of the 
solvent. Evidently, the mechanism causing the 
increased relaxation time with increasing C is 
dependent only on the number of lactoglobulin 
molecules in the solution. The value of 7 obtained 
by extrapolation to C=0 of the curve given in 
Fig. 8, is very close to 7.5+0.5X10-8 sec., in 
agreement with the value obtained by Ferry and 
Oncley.! 

It should be remarked that the absolute values 
given for the relaxation time in Fig. 8 are subject 
to considerable uncertainty due to the lack of an 
adequate theory for the local field in polar liquids. 
The results shown in Fig. 8 are based on the 
assumption that the Onsager" theory of the local 


4 L. Onsager, J. Am. Chem. Soc. 58, 1486 (1936). 
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field is valid for the B-lactoglobulin solutions. If 
on the other hand, the Clausius-Mosotti value 
for the field were assumed valid, the values of + 
would be reduced by a factor (€.+2)/(€o+2). 

The form of the dielectric dispersion is par- 
ticularly important with reference to questions 
regarding the shape of the lactoglobulin molecule. 
Ferry and Oncley! studied the dispersion of lacto- 
globulin solutions containing up to 7.7 g/l at 
temperatures of 0° and 25°C. Over the concen- 
tration range studied by them, Ferry and Oncley 
found no significant differences in the form of the 
dispersion or in the magnitude of the relaxation 
time. Comparison of their experimental disper- 
sion curves with curves calculated from the 
Debye theory for a single relaxation time (spher- 
ical molecule) and on the basis of two relaxation 
times for an spheroidal molecule showed reason- 
ably good agreement of the 25° data with either 
of the theoretical curves. For the 0° data, how- 
ever, the agreement was considered better in the 
case of the theoretical curve for the spheroidal 

. model. Shack, studying the dielectric adsorption 
of 8-lactoglobulin, obtained data indicating the 
existence of a single relaxation time.!® 

It has been shown above that the dispersion 
data obtained in the present investigation can- 
not be represented by the Debye theory for a 
single relaxation time. Moreover, it appears to be 
of doubtful significance to attempt an interpre- 
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Fic. 7. Frequency variation of dielectric constant and 
loss factor, according to Eqs. (3) and (4), for solution 
containing 94.0 grams per liter of B-lactoglobulin dissolved 
in 2.5 molar aqueous glycine at 20.0°C. 


16 J. L. Oncley in Chapter 22 of E. J. Cohn and J. T. 
Edsall, Proteins, Amino Acids, and Peptides (Reinhold 
Publishing Corporation, New York, 1943). 
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tation based on an assumed spheroidal or ellip- 
soidal molecule with two (or three) relaxation 
times. This follows from the fact that the form of 
the dispersion is not constant but depends 
markedly upon the lactoglobulin concentration 
(see Fig. 5), and if the observed dispersion is to be 
interpreted on the basis of more than one relax- 
ation time then it would appear necessary to 
assume a broad and approximately continuous 
distribution of relaxation times. Otherwise as 
Cole and Cole* have indicated, the ¢’’ versus ¢’ 
locus (see Figs. 2-4) would be required to have 
90° angles of approach to the ¢’ axis for zero and 
infinite frequencies. 

Ferry and Oncley did not report values of e’’. 
Such values are tedious to calculate, in order to 
prepare a plot of e”’ versus é’ for their data. How- 
ever, it has been shown by trial that Eq. (3) fits 
the data of Ferry and Oncley satisfactorily if @ is 
chosen equal to 0.04+0.01. It would appear 
impossible, therefore, on the basis of existing data 
for dilute solutions of lactoglobulin to distinguish 
between an interpretation based upon an assump- 
tion of two relaxation times (spheroidal model) 
and an interpretation involving some factor 
empirically described by a. However, for the 
more concentrated solutions, the assumption of 
a spheroidal or ellipsoidal molecule (two or three 
relaxation times) appears inadequate to explain 
the data. Whether §-lactoglobulin at infinite 
dilution should be represented by more than one 
relaxation time cannot be concluded from existing 
data. 

On the other hand, for concentrations of lacto- 
globulin of the order of 1 percent or more it seems 
clear that there is involved a symmetrical distri- 
bution of relaxation times about a central or 
most probable value. On the assumption that 
this interpretation is correct the actual distri- 
butions can be calculated by the use of the 
distribution function given by Cole and Cole.* 

There remains, however, the problem of ac- 
counting for the mechanism giving rise to the 
assumed distribution of relaxation times. Many 
possibilities have been suggested'*!” but it is 
difficult in view of the present state of develop- 
ment of the theory of polar solutions to do more 


* See reference 8, Eq. (13). 


16 W. A. Yager, Physics 7, 434 (1936). 
17W. Kauzmann, Rev. Mod. Phys. 14, 12 (1942). 
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Fic. 8. Variation of relaxation time, t=T0(qw/n.), with concentration of 8-lactoglobulin C. 


Data are shown for solvents containing 


than acknowledge the reasonableness of such a 
behavior. It must certainly be agreed that a 
solution of highly polar protein molecules dis- 
solved in a polar medium (aqueous glycine) offers 
many opportunities for interactions which un- 
doubtedly would be reflected in the relaxation 
time. In fact, it would be notable if they were not. 
The data presented in Fig. 5 show clearly that in 
a solution containing a fixed concentration of 
glycine, the distribution parameter a depends 
only on the number of lactoglobulin molecules in 
solution. It also shows that for a fixed concen- 
tration of lactoglobulin, a decreases as the glycine 
concentration increases. This behavior suggests 
that glycine acts to decrease the interaction 
between lactoglobulin molecules and parallels the 
well-known action of glycine to increase the 
solubility of lactoglobulin in water. 

The foregoing discussion has shown that the 
character of the dielectric constant dispersion of 
B-lactoglobulin solutions is affected by the con- 
centration of protein in the solution and by the 
concentration of the glycine in the solvent. It has 
also shown that these concentration effects can be 
formally described in terms of the parameters a 
and 79 from Eq. (2). 

As is well known, the physical properties of a 
number of proteins in solution apparently change 
as the concentration of the protein is increased 
above 0.5 to 1 percent. For example, the sedi- 
mentation constant for most proteins decreases 
with increasing concentration, markedly for 
nucleohistone, even at 0.1 percent.! On the other 
hand, lactoglobulin in inorganic salt solution 


18 R. O. Carter, J. Am, Chem. Soc. 63, 1960 (1941). 





0.48, 1.5, and 2.5 moles glycine in water. 


showed a barely significant decrease in sedi- 
mentation constant even at 4.47 percent protein.'® 
It would appear then that lactoglobulin aggrega- 
tion is not responsible for the observed effects of 
concentration on dielectric properties unless the 
behavior of this protein in glycine differs from its 
behavior in inorganic salt solution. However, a 
difference in behavior may be considered likely 
since certain proteins, readily dissolved by neutral 
salts, are not appreciably dissolved by amino 
acids such as glycine.”° Actually, a dipole such as 
glycine might be expected to be less efficient in 
suppressing interaction between §-lactoglobulin 
dipoles than an ionic salt such as sodium chloride, 
which can supply positive and negative ions to 
the negative and positive poles of the dipole, 
respectively, by ionic interchange, whereas the 
glycine carries along physically an equal and 
opposite charge. 

Ferry and Oncley”' measured the dielectric 
dispersion of serum albumin in concentrations 
ranging up to 5.1 percent and for serum pseudo- 
globulin in concentrations ranging up to 2.4 
percent. Although the dispersion curves did not 
correspond to the spherical Debye curve, no 
marked change in the shape of the dispersion 
curves or in the relaxation time occurred with 
increasing concentration. On the other hand, the 
data of Arrhenius” for hemoglobin (concentra- 
tions approximating 3 percent, except for one 


19K, O. Pedersen, Biochem. J. 30, 964 (1936). 

20 E, J. Cohn, page 619, reference 15. 

J. D. Ferry and J. L. Oncley, J. Am. Chem. Soc. 60, 
1123 (1938). 

2S. Arrhenius, Nova Acta Reg. Soc. Sci. Upsaliensis 
[IV] 12, No. 5 (1940). 
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TABLE I. Dielectric constants of lactoglobulin liquid 
crystals (32.4 percent protein). 











Frequency, v 
(megacycles) Dielectric constant, ¢ 

13.7°C 16.6°C 20.0°C 
5.00 69.0 69.3 70.9 
4.00 70.3 71.4 
3.00 71.6 69.9 70.9 
2.00 ta 71.7 
1.00 80.6 76.0 73.9 
0.800 79.6 75.0 
0.600 83.0 81.8 i, 
0.500 84.0 80.4 
0.400 84.8 85.8 83.9 
0.300 87.8 87.1 
0.200 93.3 93.5 94.2 
0.100 109.0 108.0 
0.0800 113.3 118.2 115.3 
0.0600 123.6 127.6 
0.0500 128.1 136.1 136.3 
0.0400 139.1 148.3 
0.0300 151.1 163.1 168.0 
0.0200 178.1 200.9 
0.0100 250.3 266.1 282.6 








TABLE II. Dielectric increments and dipole moments for 
B-lactoglobulin dissolved in.aqueous glycine. Dipole mo- 
ments based on Eq. (6). 














Lacto- 
Glycine globulin 
concentra- concentra- (=) 
*tion tion 
(moles) C (g/l) (€0 —€a) Cc # (e.8.u.) 
0.48 3.48 6.5 1.87+0.3! 784 X 10718 
10.4 19.2 1.84+0.1 775 
20.9 33.0 1.58+0.05 575 
LS ia 25.5 2.01+0.07 810 
iz 38.2 1.80+0.05 766 
38.2 67.3 1.76+0.03 758 
63.7 118.2 1.86+0.01 781 
2.5 11.8 25.8 2.18+0.09 844 
23.5 51.9 2.21+0.04 850 
47.0 93.0 1.98+0.02 805 
70.5 126.5 1.80+0.01 766 
94.0 167.3 1.78+0.01 764 








1 Estimated probable error in dielectric increment. 


experiment where the concentration was 9 per- 
cent) compared with the data of Oncley‘ for 
hemoglobin (concentration 0.4 to 1.4 percent) ° 
suggests a concentration effect similar to that re- 
ported here for lactoglobulin. Arrhenius reported 
a relaxation time approximately twice that re- 
ported by Oncley. In fact, such a concentration 
effect would bring their results into accord. The 
results of Johnston and Williams,” however, 
show a concentration effect for nitrobenzene dis- 


%J.H. L. Johnston and J. W. Williams, Phys. Rev. 
34, 1483 (1929). 
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solved in oil. The proportional increase in the 
relaxation time observed by these authors for a 
change in concentration from 4 to 8 percent is 
about the same as that found for lactoglobulin for 
the same change in concentration. Unfortunately, 
the shapes of the dispersion curves for nitroben- 
zene are uncertain because it was not possible to 
make measurements in the high frequency region. 

The concentration effects found for the die- 
lectric properties of lactoglobulin solutions ap- 
parently reflect the effects of concentration on the 
interaction between individual lactoglobulin 
molecules. However, since the sedimentation 
constant is so slightly affected by concentration, 
the following questions arise: (1) are the dielectric 
properties of lactoglobulin more sensitive to con- 
centration effects than other physical properties 
(sedimentation constant, diffusion constant, 
osmotic pressure, etc.), and (2) are the properties 
of lactoglobulin in glycine solution markedly 
different from the properties in organic salt 
solution ? 


Dispersion—Liquid Crystals 


The dielectric constants of B-lactoglobulin liq- 
uid crystals are presented in Table I. No values 
of the loss factor were obtained for the liquid 
crystals owing to the fact that measurements of 
thea.c. conductivity were not made at sufficiently 
low frequencies to determine Go and thus permit 
calculation of e’’ according to Eq. (1). The lack of 
values of ¢’’ precluded a test of Eq. (2) on the 
basis of the e’’ versus ¢’ plot as was done for the 
B-lactoglobulin solutions in glycine. By successive 
trials, however, Eq. (3) has been fitted to the 
data given in Table I for the liquid crystals at 
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Fic. 9. Frequency variation of dielectric constant, ac- 
cording to Eq. (3), for liquid crystals containing 324 grams 
per liter of B-lactoglobulin. The data shown are for 13.7°C. 
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Fic. 10. Dependence of infinite frequency dielectric constant €.. on lactoglobulin 
concentration C. The curve was calculated according to Eq. (8) with v=0.75 cc 


per gram. Ace. = (€s—€x). 


13.7°. The resulting curve is shown in Fig. 9. It is 
estimated that a change of +5 percent in any one 
of the four parameters €o, €~, 7) and a of Eq. (3) is 
sufficient to destroy the agreement between the 
experimental data and Eq. (3). It is noteworthy 
that the values of r>=7.6X10-* sec. and a=0.40 
for the liquid crystals are considerably larger than 
the corresponding values for the most concen- 
trated solutions of 8-lactoglobulin in glycine. 

The dielectric constants of the liquid crystals 
at 16.6° and 20.0°C differ only slightly at most 
frequencies from the values obtained at 13.7°C. 
The greatest differences occur at the lowest fre- 
quencies, as is shown in Table I, and amount to 
a maximum increase in eat 0.01 Mc of 6 percent 
(at 16.6°C) and 12 percent (at 20.0°C). 

In view of these small differences in ¢’, no 
attempt was made to fit Eq. (3) to the data for 
the higher temperatures. No claim of great accu- 
racy is made for the values of €0, €«, To, and @ 
found for the 13.7° data. However, the values 
given are consistent as to (1) the width of the 
frequency region of the dispersion of ¢’, and (2) 
the order of magnitude of the dielectric incre- 
ments and the relaxation times found for the 
more concentrated solutions of lactoglobulin in 
glycine. 


The liquid crystals furnish an example of a 
system of 8-lactoglobulin in which considerable 
molecular interaction may be presumed to occur. 
The extreme values of a and 7» found for this 
material would appear to support the suggestion 
made above that the type of dielectric dispersion 
observed for the glycine solutions of 6-lacto- 
globulin is caused by interaction between lacto- 
globulin molecules. 


Dipole Moment 


The dipole moment of the 8-lactoglobulin mole- 
cule can be estimated from the dielectric incre- 


ment (€9—€.) by means of the empirical equation, 
described by Oncley,‘ 


u=A[ M(eo—€x)/C]}. (6) 


Here yu is the dipole moment in e.s.u., A is a 
constant equal to 2.8610-'§ at 20°C, and M is 
the molecular weight (40,000 for 6-lactoglobulin). 
A summary of the values of the dielectric incre- 
ment obtained from the circle diagrams of 
Figs. 2—4 is given in Table II. Values of the dipole 
moment obtained from Eq. (6) are also shown. 
An examination of the dipole moments in 
Table II reveals an apparent tendency for the 
dipole moment to increase with increasing con- 
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centrations of glycine in the solvent and for the 
moment to decrease with increasing lactoglobulin 
concentration. As is indicated in Table II the 
dielectric increments for the lower concentrations 
of lactoglobulin are subject to relatively large 
probable errors. It was noted previously that the 
observed dielectric constants are subject to a 
probable error of about +0.5 unit, and in addi- 
tion to this possible uncertainty, the possibility 
of errors in values of €9 and e,, which were ob- 
tained by extrapolation (Figs. 2-4), must be con- 
sidered. If the value of » obtained for C= 20.9 is 
discarded as being entirely too small, the other 
values of » may be represented by a mean value 
of 791 X10-'§ e.s.u. with an average deviation of 
+26X10-8 e.s.u. This value of u is somewhat 
larger than the value reported by Ferry and 
Oncley.! These authors obtained values of u equal 
to 730 and 770 at 25 and 0°C, respectively. By 
interpolation 1.=738 at 20°C. A possible expla- 
nation for the lower value of u obtained by Ferry 
and Oncley may lie in the method of obtaining ¢€ 
used by them. Values of ¢€9 were taken as the 
average of ten values of ¢e’ in the frequency range 
0.010 to 0.080 Mc. This procedure is based upon 
the fact that they were able to detect no disper- 
sion of the dielectric constant for frequencies less 
than 0.10 mc. It appears that any error intro- 
duced by their procedure would tend to give too 
small values of ¢9 and thus of yu. 

The values of €. used by Ferry and Oncley were 
calculated on the assumption that the volume of 
the solution occupied by the lactoglobulin mole- 
cules does not contribute appreciably to the 
dielectric polarization at frequencies lying beyond 
the dispersion region. As will be shown below, 
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values of €, obtained in this manner agree well 
with the values of €. obtained by extrapolation 
using the circle diagram method. 

If it is assumed that the high frequency die- 
lectric constant e, of a protein solution is con- 
tributed principally by the solvent molecules and 
that the volume occupied by the protein has a 
high frequency dielectric constant of unity, then 


(€x —€s) = (es —1)vC/1000 (7) 
or 
(€~ —€s)/(€s— 1) =vC/1000, (8) 


where v is the specific volume of the protein and 
e, is the dielectric constant of the solvent. This 
equation is subject to question from the lack of 
the true high frequency dielectric constant of the 
protein and of the degree of hydration of the pro- 
tein in solution. The equation has been fully 
discussed by Oncley.* A graph of Eq. (8) is given 
in Fig. 10. The graph has been drawn for 
v=0.75 cc/g. On the graph are plotted values of 
(€2—€s)/(€s—1) for the solutions of B-iactoglobu- 
lin and for the 6-lactoglobulin liquid crystals. The 
extrapolated values of ¢. for the glycine solutions 
were obtained from the high frequency extrapola- 
tion of the circle diagrams of Figs. 2-4. The value 
of «. for the liquid crystals was obtained from 
Fig. 9 as discussed previously. 

It is seen that within the estimated limits of 
error the values of ¢€. obtained by extrapolation 
by the circle diagram method and by means of 
Eq. (8) are in good agreement. 


* J. L. Oncley, in Chapter 22, of J. T. Edsall and E. J. 
Cohn, Proteins, Amino Acids, and Peptides (Reinhold 
Publishing Corporation, New York, 1943). 








1e 


of 
on 
of 


old 








THE JOURNAL OF CHEMICAL PHYSICS 


VOLUME 





12, NUMBER 11 NOVEMBER, 1944 


Note on Becker’s Theory of the Shock Front 


L. H. THomas* 
Ballistic Research Laboratory, Aberdeen Proving Ground, Maryland 


(Received October 16, 1944) 


Becker, in a classical paper, concludes that weak shocks are many free paths thick and may 
be treated by ordinary hydrodynamics including the effects of viscosity and thermal con- 
ductivity of the medium; that moderately strong shocks have a thickness of the order of a 
free path and must be treated by a direct attack on the relevant Boltzmann equation; while 
violent shocks have a thickness small compared with a free path and even the Boltzmann 
equation is no longer applicable. It is shown here that this last conclusion rests on an oversight; 
the thickness of a shock front is always at least of the order of magnitude of a free path, and 
it is to be expected that the Boltzmann equation can be applied even for the most violent shocks. 





INTRODUCTION 


ISCOSITY and heat conduction are often 

neglected in the equations of motion of a 
compressible fluid. The researches of Earnshaw! 
and Riemann? showed that a continuous solution 
of the resulting equations may be impossible for 
longer than a finite time if an initial disturbance 
moves from high to low pressure. The velocity 
of sound being greater at greater temperature 
produced by adiabatic compression, a disturb- 
ance on the high pressure side moves faster than 
on the low, and the disturbance becomes more 
and more concentrated, until finally no con- 
tinuous solution is possible. Physically, a shock 
wave has been built up. The subsequent motion 
of the shock front may be obtained by regarding 
it as a surface of discontinuity across which 
there is continuous flow of mass and momentum, 
leading to equations given by Stokes’ and later 
by Rankine.‘ In order that energy also should be 
conserved the entropy must increase across the 
front leading to a further condition owing to 
Huguoniot.’ (A rarefaction shock, across which 
entropy would have to diminish, is therefore not 
possible.) 

It was explained by Rankine and later by 
Rayleigh® and by Taylor’ that if viscosity and 
heat conduction are taken into account, the 
shock wave would never become an infinitely 


* On leave of absence from Ohio State University. 
'S. Earnshaw, Phil. Trans. 150, 133 (1858). 

2 B. Riemann, Gott. Abh. 8, 43 (1858). 

°E. E. Stokes, Phil. Mag. 3, 23, 349 (1848). 

*W. V. M. Rankine, Phil. Trans. 160, 277 (1870). 
* A. Hugoniot, J. Ecole Poly. 58, (1889). 

* Lord Rayleigh, Proc. Roy. Soc. A84, 371 (1910). 
7G. 1. Taylor, Proc. Roy. Soc. A84, 371 (1910). 
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sharp discontinuity, but would have a finite 
thickness (3X10-° cm for moderately strong 
shock waves in air). Becker,* in a classical paper, 
obtained the exact solution of these equations 
for a one-dimensional stationary shock wave in 
a simple gas. He found that for a violent shock 
(p2/p1=8, where pz and p; are the pressures on 
opposite sides of the shock), the thickness for a 
gas was of the order of a mean free path (10-5 
cm). The equations of hydrodynamics may be 
derived from kinetic theory only provided that 
large macroscopic changes do not occur over a 
distance of the order of a free path. Becker 
therefore concludes that these equations would 
not apply within a moderately violent shock, 
and that recourse must then be had to a direct 
solution of the Boltzmann equation of kinetic 
theory. 

Becker finds further that for a very violent 
shock (p2/p1= 2000) the thickness of the shock 
wave on the basis of these (inapplicable) equa- 
tions is so much smaller yet as to approach the 
average intermolecular distance (310-7 cm) 
and he inferred that the same order of magnitude 
holds for liquids. He therefore doubts even the 
applicability of the Boltzmann equation to very 
violent shocks. 

Unfortunately Becker overlooks the increase 
of the coefficients of thermal conductivity and of 
viscosity with increasing temperature and pres- 
sure, using the values for normal temperature 
and pressure as constants. If this oversight is 
corrected the hydrodynamical equations never 
lead to a thickness small compared to the length 


8 R. Becker, Zeits. f. Physik 8, 321 (1922). 
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of a free path. Thus the ordinary approximation 
of the kinetic theory of non-uniform gases is 
valid for the less violent shocks and leads to the 
equations of hydrodynamics without viscosity 
and heat conduction as a first approximation, in- 
cluding these in the second approximation, in 
which approximation Becker’s solution is ob- 
tained. The coefficients of the terms required for 
the third approximation have been given by 
Burnett.® The inclusion of these terms leads to 
an even greater thickness, but the numerical 
values indicate that the series of approximations 
may not converge for very violent shocks, for 
which the departure from a Maxwellian dis- 
tribution of velocities at points within the shock 
front may be large. However, methods similar 
to those of Milne!® and Hopf" for the transfer of 
radiation through the solar photosphere are still 
available for solving the Boltzmann equation. 


THE EQUATIONS OF HYDRODYNAMICS AS 
OBTAINED FROM KINETIC THEORY” 


For steady motion in one dimension (x) the 
equations of continuity, momentum, and energy, 
are obtained from kinetic theory in the form 


dp du 
— 7%, (1) 
d du d 
——- ay (2) 
d du d d 
—* a (3) 


where u is the mean molecular velocity in the 
x direction, p is the density, E is the energy per 
unit volume referred to axes moving with 
velocity u, p is the mean rate of flow of momen- 
tum in the x direction per unit area perpendicular 
to the x direction referred to these axes, and q 
is the mean rate of flow of energy defined 
similarly. 

®E. S. Burnett, Proc. Lond. Math. Soc. 39, 385 (1935). 

10F, A. Milne, M.N.R.A. S. 88, 493 (1928). 

il E. . Hopf, “Mathematical problems of radiative equilib- 
rium,’’ Cambridge Tracts in Mathematics and Mathematical 
rt + (Cambridge University Press, Cambridge, England, 

12S. Chapman and T. E. Cowling, Mathematical Theory 


of Non-Uniform Gases (Cambridge University Press, Cam- 
bridge, England, 1939), Chapters 3, 7, 15. 


These equations have the first integrals® 


pu=a, (4) 
au+p=b, (5) 
uE+ub—tau?+q=c, (6) 


expressing that the rates of flow of mass, mo- 
mentum, and energy, in the x direction, are the 
same for all values of x. 

These integrated equations must hold gener- 
ally and show that, provided there exists a solu- 
tion of the Boltzmann equation corresponding 
to steady flow from uniform steady motion for 
large negative x to uniform steady conditions for 
large positive x, the relation between these 
uniform conditions is that given by the Rankine- 
Hugoniot equations to which these equations 
lead. 

The usual first approximation of kinetic theory 
gives, for a gas with energy only of molecular 
translation, 


gmkT, p=p®=nkT, q=0, (7) 


where 7 is the number of molecules per unit 
volume, T is the absolute temperature, k is 
Boltzmann's constant, and m is the mass of a 
molecule. 

The methods of Chapman and Enskog, using 
p=nm, E=3nkT (the latter essentially as a 
definition of T), lead to a solution of Boltzmann’s 
equation by successive approximation in series 
giving 


p=nm, E= 


p=pO+p%+4+pO4 rey 





(8) 
g=qU+qO+---, 
where, in the one-dimensional case,” 
O=nkT, po a ‘sll (9) 
gr=aai, pre--—2—, ("= -Ar— 
3 dx dx 
uw? du2du iT 2d s1 dp 
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Pp dx3dx pl 3dx\pdx 
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uw? dudT w(2d/ du 
oneal S(t 
pl dx dx pl i3dx\ dx 


uw dp 2 du 


op dx 3dx 

pw? d s2du uw? dT 2du 
+o—-—(“—) 49 

p dx\3 dx 


———-—, (10) 
pI dx 3 dx 
where uz is the coefficient of viscosity, \ the coef- 
ficient of thermal conductivity, the forms of 
which have been computed for various laws of 
interaction between the molecules. 
The coefficients #1, ---, Js, and @, +++, @, 
have been evaluated for molecules attracting 
according to an inverse fifth power law :° 








31=75/8, .=45/8, 33=—3, 

34=3, 3:=6, (11) 
&1=10/3, a2=2, a=3, a =0, 

&e=3, Ge=8. (12) 


The coefficients @1, ---, @s, have also been 
evaluated for rigid elastic spheres. 

Successive approximations to the solution of 
the Boltzmann equation may now be ob- 
tained from Eqs. (4)-(6), by using only the 
terms p-+p and g®; or only the terms 
pO+p%M+p@ and g+q®; and so on; and 
subject to the convergence of this process, a 
solution of the Boltzmann equation of the desired 
form would be obtained. . 

For a gas made up of molecules with internal 
degrees of freedom this is not entirely correct; 
energy does not partition itself immediately into 
the internal degrees of freedom. Bethe and 
Teller’ have given arguments which show that 
while energy passes rapidly to rotation of the 
molecules, it passes only slowly to internal vibra- 
tion. Through the shock front itself, when of 
only a few mean-free-paths thickness, we need 
take account only of the molecular rotation. 
The shock front is then followed by a region in 
which the energy is gradually distributed until 
the equilibrium state is reached. , 

Through the shock front it would still be 
necessary, strictly, to take account of the delay 


“H. A. Bethe and E. Teller, ‘‘Deviations from thermal 
equilibrium in shock waves’ (unpublished). 


in exchanging with molecular rotation. In the 
only such case for which the theory of non- 
uniform gases has been worked out, that of 
rough spherical molecules, the first approxima- 
tion departs little from what would be obtained 
by taking the equations (4)—(6), with 


kN 
E=4NnkT, \= fu, (13) 


m 


where WN is the number of effective degrees of 
freedom (3 for monatomic, 5 for diatomic, and 
6 for polyatomic molecules) and f is a constant 
(~}(9y—5) where y=(N+2)/N is the ratio of 
specific heats. 

Finally, when we require the variation of u 
with temperature, we have, for elastic spheres 


of diameter o, 
v (kmT)} 
eli” ’ (14) 


ri = g? 





where v=0.998, while the mean free path / is 
given by 
l=1/v2rno?. (15) 


The variation of » with temperature for actual 
gases is given more nearly by a formula of Suther- 
land’s type, but this gives variation of the above 
form at high temperatures. 


BECKER’S THEORY OF THE THICKNESS OF 
SHOCK FRONTS® 


If in Eqs. (4)—(6), we replace p by p™+p™, 
i.e., nkT—(4/3)y(du/dx), and q by q®, i.e., 
—\(dT/dx), they may be regarded as equations 
to determine du/dx and dT/dx in terms of u 
and 7. Making the substitutions 


b m b? 
u=-w, T=——0 (16) 
a k a 
and writing 
. kN 
a 


these become Becker’s equations 


14 dw 0 
~<a |, 
a3 dx w 
(17) 


1 dé 
—fu—=0—6[(1—w)?+a], 
a dx 





where 
1 2ac 
a=—-l, 
N b? 
which give 


do 4 #4 6L(1—w)? +a ] 
= (18) 
de 3f- w*—w+0 





to determine the variation of @ with w across the 
shock. 

Following Becker the solution of this equation 
for which x extends from —* to ~& extends 
from 6;5=w1(1—a) to 62=we(1 —ws), where 


N+2 2ac/b? 
wi to2=—, wwe= . 
N+1 





In case 
“ 2 
— 1 +—, 
+ N 
this special solution has the simple form 
1 2ac 
saat —-u') 
~ N42 
O3wW9? — Oow 1" O.— A; 
= _ w. (19) 


wwe? wy — a2? 





For air N=5, f=1.95 (observed); if we take 
f=1.86, following Becker, we have this simple 
case, and we may expect this change of f to make 
only small differences in the result. 

Thus far, according to Becker: To find the 
thickness of the shock front, however, for violent 
shocks, for which there is considerable change of 
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Fic. 1. Definition of ‘‘thickness’’ of the shock front. 


pressure and temperature, it is essential to take 
account of the increase of viscosity and thermal 
conductivity with temperature. 

Using (14) and (15), if we measure ‘“‘distance’’ 
in terms of ‘‘free paths,”’ i.e., in terms of s where 


ds 
—=!l, (20) 


we replace Becker’s equations (17) by 


47/2\3 /0 dw 0 
-(- v— —=w+--l, 
3\r w ds w 
(21) 
2\2 6 d0 
or: 
T w ds 


“fo 9] 


We may measure the thickness of the shock 
front by the whole change in w divided by the 
greatest slope dw/ds (Fig. 1). 





ds 472\3 ((w1 
oxo (=) 35)" 


This thickness, as well as the maximum value 
of p/p in the front, 


= 4p “s 
2  —- 3 pdx} max 


| (w? — w) (wi + we) 


1 23 
w?(1—wi— a ( ) 


we) + W1W2 max 





is given for shocks of various intensity in 
Table I. 


w1+we— 1)w?) |} 
(wit ) = . (22) 


(w1 —w) (w— we) 





If now the expressions (17) for dw/dx and d6/dx 
are modified by including p® and g® evaluated 
from the above solution on the right-hand sides, 
they give 


14 dw p 
SE (o)/(HB) 
a3 dx pe 


“ice=(0-s1-w)-t0t3/(14%2), 


(1) 


(24) 
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TABLE I. 

p2 pe Te pm 

_— - —_ Thickness 

pi pi T1 (free paths) po 
1.0000 1.0000 1.0000 0 0.0000 
2.6661 ~ 4.4981 1.6871 3.98 0.3460 
3.7899 9.8365 2.5954 3.08 0.7198 
4.6382 19.700; 4.2474 2.25 1.0967 
5.2931 43.3876 8.1998 1.98 1.4854 
6.0000 0 2% 1.74 2.5000 








where, from Burnett’s formulae, (9)—(12), 


p/p = p/p, 
and , 
g/g =10p/p., 


We see that, while there is no guarantee of con- 
vergence, the general effect of the third-order 
terms is to make the front thicker. 

Lastly, in the notation of kinetic theory, the 
Boltzmann equation, the solutions of which we 
are approximating, is 


aff ff rn -s10sbabdeiosn, 
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to determine f(x,u,v) and, if the method of 
Chapman and Enskog fails to converge, this may 
be rewritten 


ssf fff saabdciosn 
=f fff riverardedo,au, 


and may be approximated by iteration using an 
already determined approximation both on the 
right-hand side and in the integral on the left, 
in the manner used in astrophysics for the 
transport of radiation in a stellar photosphere.!° !! 


CONCLUSIONS 


All shock waves in air but the weakest are a 
few free paths thick. If the shock is strong the 
approximation process of the theory of non- 
uniform gases probably does not converge, 
because the distribution function of molecular 
velocities within the front departs too far from 
the Maxwellian. There is, however, no reason to 
doubt that the Boltzmann equation possesses a 
suitable solution. 
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The Charge Effect in Photographic Development: Role of the Gelatin* 


T. H. JAMEs 
Eastman Kodak Company, Rochester, New York 


(Received October 2, 1944) 


The approach of negatively charged developer ions to the surface of a silver bromide grain 
is opposed by charge barriers originating both with the gelatin and with the adsorbed bromide 
ions. It is shown that the effect of neutral salt upon the rate of development by negative ions 
results primarily from a depression of the gelatin charge barrier. The effect of salt upon the 
bromide barrier is relatively unimportant. On the other hand, a change in the bromide barrier 
is mainly responsible for certain other kinetic effects, such as the induction period in develop- 
ment and the accelerating action of phenosafranin. 


COLLOIDAL particle of silver bromide, 

suspended in a solution containing excess 
soluble bromide, is negatively charged by the 
adsorption of excess bromide ions to the surface. 
This charge acts to concentrate positive ions in 
the solution near the surface and to oppose the 
approach of other negative ions to the surface. 


*Communication No. 994 from the Kodak Research 
Laboratories, 


The resulting electric double layer consists on 
the one hand of the adsorbed bromide-ion layer, 
together with a more or less firmly held water 
shell, and on the other hand a diffuse ionic 
atmosphere in which positive ions predominate 
near the solid surface but decrease in pre- 
dominance as the distance from the surface 
increases. 

If the system is placed in an electric field, the 
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particle moves through the solution. The mobility 
of the particle is determined by the potential at 
the surface of shear between the particle and 
the bulk solution, i.e., the zeta-potential. For a 
given concentration of bromide ions in the 
adsorbed layer, the potential at any point in the 
solution near the adsorbed layer depends 
markedly upon the ion concentration in the 
solution. The higher the ion concentration, the 
greater will be the crowding of positive ions near 
the solid surface, and hence the lower the poten- 
tial in this region or at the shear surface itself. 

In the presence of gelatin, the silver bromide 
particle becomes coated with an adsorbed layer 
of gelatin which is tenaciously held to the sur- 
face.! It is not known to what extent the gelatin 
displaces the adsorbed bromide and water layer, 
and to what extent it simply forms an additional 
outer coating. However, the behavior of the 
particle in an electric field is now determined by 
the gelatin and not the bromide layer. On the 
positive side of the gelatin isoelectric point, the 
particle moves as a positively charged body; on 
the negative side, it moves as a negatively 
charged body.? The surface of shear is between 
the gelatin layer and the solution. 

The sensitive layer of the normal photographic 
material consists fundamentally of tiny crystals 
of silver halide (the grains) embedded in a layer 
of gelatin. These crystals are prepared and aged 
in the presence of excess halide and gelatin, so 

1S. E. Sheppard, R. H. Lambert, and R. L. Keenan, 
J. Phys. Chem. 36, 174 (1932). 


2 Cf. L. F. Tice and W. G. Batt, Am. J. Pharmacy 109, 
29 (1937). 
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that adsorbed layers of each should form as in 
the case of the colloidal silver bromide particle 
just discussed. The adsorbed gelatin layer is now 
surrounded, however, by bulk gelatin instead of 
by solution. 

In normal “‘chemical,”’ or contact development 
of the photographic layer, the reaction takes 
place at the interface between the silver halide 
and silver (either from latent image or previous 
development). The developing agent accordingly 
must penetrate to the grain surface before reac- 
tion can take place. Many of the developing 
agents are active only in the form of negatively 
charged ions. In such cases, it has already been 
shown’ that the potential in the vicinity of the 
grain surface exerts a pronounced influence upon 
the kinetics of development. 

In previous investigations of this charge effect, 
development was always carried out in alkaline 
solution at a pH well above the isoelectric point 
of the gelatin. Accordingly, both the bromide 
layer and the gelatin could contribute to a charge 
barrier resisting the approach of the developer 
ions to the surface of the silver bromide grain. 
The existing data on the effect consisted largely 
of measurements of the change in rates of devel- 
opment produced by the addition of bromide, 
various dyes, and neutral salts to the developer 
solution. Certain other development charac- 
teristics, particularly the existence and mag- 
nitude of the induction period, were also shown 
to depend upon the charge of the developer ion 
rather than upon the chemical nature of the 
agent. These data seemed to be largely ex- 
plicable on the basis of the bromide charge 
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Fic. 2. Salt effect in terms of ratio of rates in presence 
and absence of 0.33 M KNO3. Curve A is for film with a 
gelatin of isoelectric point=4.9; curve B is for film with 
a gelatin of isoelectric point =8. 


3T.H. James, J. Phys. Chem. 43, 701 (1939). 
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alone,** but a hint of some action of the gelatin 
was contained in observations on the behavior of 
the ferro-oxalate developer in acid solution.’ 
Recently, however, it has been shown® that 
negatively charged gelatin layers exert a definite 
restraining action on the penetration of nega- 
tively charged dye ions, and hence might be 
expected to exert a similar action upon the 
penetration of developer ions. The resistance to 
the penetration of the dye ions was largely 
removed by the addition of large amounts of 
neutral salt, and Sheppard has suggested® that 
a similar action can account, at least in part, for 
the increase in rate of development produced by 
adding neutral salts to solutions of ionic develop- 
ing agents.’ 

In the present paper, it will be shown that a 
change in the electric barrier of the gelatin is 
almost entirely responsible for the acceleration 
produced by neutral salt. The slight effect of salt 
upon the bromide barrier can be measured by 
working at the isoelectric point of the gelatin. 
By following this procedure further, it can be 
shown that the bromide barrier is primarily 
responsible for the induction period and for 
certain other kinetic effects, such as the acceler- 
ating action of phenosafranin. 


MATERIALS AND PROCEDURE 


Two distinct photographic materials were em- 
ployed. The first (A) was a normal motion- 
picture positive film prepared from gelatin which 
had an isoelectric point at about 4.9. The second 


(B) was a positive-type film prepared from an . 


acid-process gelatin which had its isoelectric 
point at about 8. The light-sensitive material in 
each case was silver bromide. 

As a developer, the bivalent ferro-oxalate is 
well suited to the purpose of the investigation. 
The oxidation product produces no complicating 


*F. Urbach, [X* Congrés International de Photographie 
(Paris, 1935), p. 432. J. E. de Langhe, Zeits. f. wiss. Phot. 
35, 201 (1936). A. Lottermoser and R. Steudel, Kolloid 
Zeits, 82, 319; 83, 37 (1938). T. H. James, J. Phot. Soc. 
Am. 9, 62 (1943). 

°S. E. Sheppard, R. C. Houck, and C. Dittmar, J. Phys. 
Chem. 46, 158 (1942). 

*S. E. Sheppard, “‘Colloid chemical aspects of photo- 
graphic development,” in Colloid Chemistry, edited by 
J. Alexander (Reinhold Publishing Corporation, New 
York, 1944), p. 472. 

7H. Liippo-Cramer, Kolloid Zeits. 15, 164 (1914). 
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Fic. 3. Salt effect on film with gelatin isoelectric point 
of approximately 8. —-O—O=pH=4.7; log E=1.45; 
—X— X =pH=12.7; log E=1.75. 


kinetic effects under the experimental conditions 
employed. The developer can be used in acid, 
neutral, and slightly alkaline solution. Thus, de- 
velopment can be carried out on either side of 


‘the isoelectric point of the gelatin. The bivalent 


iso-ascorbic acid ion was used as developer in 
the highly alkaline region, since in this case also 
the oxidation product does not complicate the 
development kinetics. *® 

The ferro-oxalate developer contained 0.005 
ferro-oxalate and 0.04 M excess potassium oxa- 
late. The pH of the solution was adjusted by 
adding the required amount of sulfuric acid for 
pH values of less than 7.5, and borax for higher 
values. The concentration of iso-ascorbic acid 
was 0.005 M in a solution containing 0.04 M@ 
potassium hydroxide. The developers were pre- 
pared from oxygen-free component solutions, and 
all development operations were carried out in 
an atmosphere of tank nitrogen which had been 
freed of oxygen. The apparatus and procedure 
were essentially those previously described.’ 
Development rates were determined in terms of 
1/t and R, where ¢ is the time required to obtain 
a fixed, small optical density at a given ex- 
posure, and R is the slope of the curve obtained 
by plotting density against time of development. 
All values are corrected for fog. Unless otherwise 
specified, the data recorded here apply to an 
exposure of log E=1.45, where E is expressed in 
meter-candle-seconds, ¢ is measured at a density 
of 0.20, and R at a density of 0.80. Similar 
results were obtained for other exposure values 
in the normal exposure range. 


8’ T. H. James, J. Am. Chem. Soc. 66, 91 (1944). 
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TABLE I. Salt effect in development at the isoelectric point. 
Ferro-oxalate developer, pH = 4.9, film A 











KBr KNOs i/t R 
0.0 0.0 0.099 0.0485 
33 (098 048 
‘67 097 046 
00067 0 034 036 
67 036 041 
0020 0 022 029 ~ 
67 025 033 








RESULTS AND DISCUSSION 


The rate of development of film A in ferro- 
oxalate decreased with increasing pH, as in- 
dicated by the solid curve of Fig. 1. The addition 
of neutral salt changed the rate of development 
in most cases, but produced no change at the pH 
corresponding to the isoelectric point of the 
gelatin (4.9). The broken curves in Fig. 1 illustrate 
the effect of potassium nitrate upon the develop- 
ment rate. On the positive side of the isoelectric 
point, the salt produces a decrease in the rate of 
development; on the negative side, the salt 
produces the more conventional increase. The 
three curves, corresponding to no added salt, 
0.33 M potassium nitrate, and 0.67 M potassium 
nitrate, cross at the isoelectric point. Quite 
similar results were obtained when R was 
plotted instead of 1/t. 

The absolute rate of development also can be 
related to the charge of the gelatin. In the 
absence of salt, increase in the positive charge of 
the gelatin increases the rate of development, and 
increase in the negative charge decreases the 
rate of development. At high salt concentrations, 
the charge effect of the gelatin is practically 
eliminated,‘ and the rate of development becomes 
almost independent of the pH. 

As a check on the results obtained with film 
A, the effect of salt on development of film B 
was measured. Figure 2 illustrates the results 
obtained. In this figure, the ratio of development 
rate in the presence of 0.33 M potassium nitrate 
to that in the absence of added salt is plotted 
against pH (curve B). The ratio is unity(no salt 
effect) at a pH of about 8.1. This corresponds 
well with the isoelectric point of the gelatin used 
in the preparation of this film. For comparison, 
results obtained with film A are plotted in the 


JAMES 


same way (curve A). The point on this curve 
corresponding to pH=12.7 is for iso-ascorbic 
acid instead of ferro-oxalate, but the extension 
of the curve follows the general shape indicated 
by the ferro-oxalate points. Thus, the specific 
properties of the individual developing agents 
do not appear to be involved in determining the 
magnitude of the salt effect in these two cases. 
Only the charge is important. 

A series of measurements of the salt effect at 
various salt concentrations was made at pH =4.7 
with ferro-oxalate and 12.7 with iso-ascorbic 
acid, using film B as the photographic material. 
(Data for film A at pH=12.7 have been given.®) 
The determinations with the ferro-oxalate de- 
veloper were made in the absence of added 
bromide, but addition of 0.00667 M potassium 
bromide to the iso-ascorbic acid developer was 
necessary to eliminate excessive fog. Results of 
the determination at the two pH values are 
given in Fig. 3 as the rate of development, 1/1, 
plotted against the square root of the potassium- 
ion concentration. The salt effect is large in each 
case, but is opposite in direction. 

All of the preceding experiments with ferro- 
oxalate were carried out in the absence of added 
bromide. Only the bromide in excess in the 
original film was present. Under these conditions, 
the evidence shows clearly that the salt effect is 
acting exclusively upon the gelatin charge, within 
the limits of experimental error. This does not, 
however, preclude an effect upon the bromide 
when the latter is in considerable excess. Ac- 
cordingly, measurements in the presence of 
added bromide were made at the isoelectric point 
of the gelatin, where the electric barrier of the 
latter is eliminated. Some results are given in 
Table I. At the higher bromide concentrations, a 


TABLE II. Effect of bromide on development by ferro- 
oxalate. Film A—no added salt. 











pH KBr 1/t R Ri 
4.76 0.0 0.118 0.054 0.45 
.000167 071 052 0.73 
.000667 044 .049 1.11 
.00167 .030 .039 1.30 
4.06 0 192 .080 0.42 
.000167 167 .087 0.52 
.000667 114 .084 0.76 
.00167 .080 075 0.94 
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THE CHARGE EFFECT 


TABLE III. Effect of phenosafranin. Film A; pH =4.06; 
KBr, 0.001667 M. 











Phenosafranin 1/t R 
0.0 0.080 0.075 
.00002 -105 .084 
.00005 .132 .099 
.00010 154 .095 
-00020 .192 .100 
.00050 455 .126 








small increase in development rate is obtained by 
increasing the salt concentration. However, the 
effect is very small in comparison with that ob- 
served at similar salt concentrations, but in 
alkaline solution well to the negative side of the 
isoelectric point. The neutral salt acceleration of 
development, therefore, is caused primarily by an 
action upon the gelatin charge barrier, even when 
considerable excess bromide is present. 

It was shown previously* that the induction 
period in development is caused by a charge 
effect. The original explanation of this effect 
associated it with a change in the adsorbed 
bromide layer as development progressed. This 
explanation remains in agreement with the 
evidence. The induction period persists even on 
the positive side of the isoelectric point. More- 
over, its characteristic sensitivity to the bromide- 
ion concentration® ® is maintained on the positive 
side of the gelatin isoelectric point. This is indi- 
cated by the data given in Table I]. It will be 
noticed that the 1/t rates, which are measured 
in the induction period region, decrease much 
more rapidly with increasing bromide concen- 
tration than do the R rates, which apply to the 
region of development beyond the induction 
period. 


®°T. H. James, J. Phys. Chem. 44, 42 (1940). 
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The relative change in 1/t with respect to R 
is indicated by the figures in the last column. 
The increase in Rt with increasing bromide is 
not as great at PH=4.06 as at 4.76, probably 
because of the tendency of the positive gelatin 
layer to counteract to some extent the kinetic 
effects of the negative bromide barrier. Similarly, 
the absolute decrease in rate of development with 
increasing bromide is not as great at the lower 
pH—a fact which can be explained in the same 
way. However, the induction period persists at 
pH=4.06, and the relative extent of the in- 
duction period still increases with increasing 
bromide-ion concentration. 

The action of phenosafranin supplies another 
indication that the bromide barrier is still 
actively retarding development well to the 
positive side of the gelatin isoelectric point. 
This substance increased the rate of development 
by all of the doubly-charged agents previously 
studied, and its action is definitely associated 
with the charge of the developer.* In the present 
investigation, it was observed that phenosafranin 
accelerates development at the isoelectric point 
of the gelatin, and indeed well to the positive 
side of this point. The data in Table II illustrate 
the effect at pH=4.06, and in the presence of 
0.001667 M potassium bromide. The acceleration 
produced by the dye is marked, and comparison 
with Table II shows that 0.0002 M phenosafranin 
almost completely erased the retarding action 
which the bromide exerted on development in 
the absence of the dye. Thus, under the present 
conditions, approximately one molecule of pheno- 
safranin offsets the retarding effect of eight 
bromide ions. Approximately the same relation 
held over the sixteen-fold range of exposures for 
which data were tabulated. 
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The Heat Capacity of Gaseous Butadiene-1,3* 
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The heat capacity of gaseous butadiene-1,3 at 300°K and at 363°K has been determined by 
the measurement of the velocity of high frequency sound in the gas over the pressure range 
from one-third to one atmos. The heat capacity of the gas at zero pressure, evaluated by 
extrapolation to zero pressure, is in essential agreement with the value calculated from experi- 
mental state data. The experimental results are compared with values calculated from spec- 
troscopic and molecular data and a marked disagreement is observed. 





INTRODUCTION 


HE present paper describes some heat 

capacity data at 300°K and 363°K ob- 
tained by measurements of the velocity of ultra- 
sound in gaseous butadiene-1,3. The purpose of 
the measurements was to obtain experimental 
gaseous heat capacity data for comparison with 
those calculated from spectroscopic and mo- 
lecular data. The agreement of C,° values 
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Fic. 1. The velocity of ultrasound in butadiene-1,3. 
@ and scale—t=25.82°C. © and right scale— 
t=90.13°C. , 


derived from the velocity of ultrasound in the 
gas at one temperature for several pressures 
constitutes a verification of the state data used 
in the computation. In this work we have used 
state data! obtained at the National Bureau of 

* The investigation reported in this paper is part of a 
joint program with the Low Temperature Section, National 
Bureau of Standards, on the C, hydrocarbons. 

1C. H. Myers, R. B. Scott, F. G. Brickwedde, and R. D. 
Rands, Jr., “Thermal properties of vapor and liquid 


butadiene-1,3,"" submitted to the Office of the Rubber 
Director (June 30, 1943). 
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Standards. The results confirm the validity of 
these data. The comparison of the experimental 
heat capacities at zero pressure with those from 
spectroscopic and molecular data will be dis- 
cussed presently. 


MATERIALS 


The butadiene used was Dow Chemical Com- 
pany material obtained from the Petroleum 
Refining Laboratory of The Pennsylvania State 
College, through the courtesy of Professor M. R. 
Fenske. 

The butadiene was dried by passing it twice 
over freshly activated alumina at the rate of 
about 1 mole per hour. The alumina was pre- 
viously baked at 300-330°C for 12 hours under 
high vacuum before each drying. The dried 
butadiene was distilled through the low tem- 
perature column of the Cryogenic Laboratory 
(12 plates). The cut used for the supersonic 
measurements had a boiling range of 0.05°K 
(normal boiling point, 268.81°K) which was 
assumed to indicate satisfactory purity. 


EXPERIMENTAL METHOD AND RESULTS 


The apparatus and method used were identical 
with those of Telfair? and the calculation was also 
based on the same equations. 

The heat capacity of the gas at infinite dilution 
was computed from the velocity of sound in the 
gas with the equation:? 

RT QR 
y-—|P+—|, 
M C, 


(1) 


where P and Q are dimensionless functions of 
2D. Telfair, J. Chem. Phys. 10, 167 (1942). 
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the virial coefficients and their derivatives, 
defined by 


pu=RT+B/v+8/0+e/e+---, (2) 
P=1+28/RTv+38/(RTv)?+--°, (3) 


1/96 1/06 : 
o-[1+_(—) + (—) tee], (4) 
Rv\ dT/, Rv\dT/, 


V is the velocity of sound, C, is the heat capacity 

of gas at constant volume under the given con- 

ditions, p is the pressure, and v the molal volume. 
Solution for C, yields 


OR 
PM 


RT 





(5) 


v 


Before beginning measurements, the gas 
chamber was baked and evacuated continuously 
for 48 hours: 36 hours at 225-250°C and 12 
hours’ at 265-290°C. This precaution was to 
insure that all volatile material was completely 
removed from the system. The initial leak rate 
on the entire apparatus was 2X10 mm/hr. 
and the system was evacuated to 5X10-* mm 
or better before beginning each series of measure- 
ments. 

In order to correct for the small temperature 
differences between points and to reduce the 
velocity measurements of each series to a 
common temperature, the differential relation 


2AV/V=AT/T 


was used. The velocity measurements in each 
series in general were taken in the order of in- 
creasing pressure. Additional sample was intro- 
duced into the chamber at the conclusion of each 
determination. 

The results are summarized in Tables | and II. 
In the tables, column 1 gives the order number 
of the point; column 2 is the pressure in cm of 
mercury; column 3 gives the experimental sound 
velocity in m/sec.; column 4 gives the tem- 
perature in °C; column 5 gives the velocity cal- 
culated for the common temperature of the 
series; and column 6 is the square of this 
quantity. 

The square of the velocity is plotted against 
pressure in Fig. 1, for determinations at 300°K 
and 363°K. It is found that a linear relationship 


TABLE I. Velocity of ultrasound (approximate frequency 
445 kc) in butadiene-1,3. 








Pressure V v2 
at 25.82°C at 25.82°C 





cm of Vexp t 
Point Hg m/sec. a m/sec. (m/sec.)? 
1A 37.2 223.90 25.98 223.84 50,104 
1B 22.6 224.85 25.75 224.85 50,567 
1C 63.0 222.16 25.68 222.21 49,377 
1D 105.5 219.56 25.88 219.53 48,193 


(AV =0.374T) 
When extrapolated to zero pressure: 


V?=51,210+30; V=226.30+0.06 m/sec.; 
C,°=17.3340.10 cal./deg./mole at 25.82°C. 








TABLE II. Velocity of ultrasound (approximate frequency 
445 kc) in butadiene-1,3. 











Pressure V 4 

cm of Vexp t at 90.13°C at 90.13°C 

Point Hg m/sec. i m/sec. (m/sec.)? 
2A 23.5 247.09 91.51 246.69 60,855 
2B 36.1 245.96 90.53 245.84 60,437 
2C 59.0 243.86 90.00 243.90 59,487 
3A 15.5 246.60 89.60 246.75 60,886 
3B 39.8 245.42 89.60 245.57 60,305 
3C 76.0 244.01 88.80 244.40 59,731 
4A 30.1 246.08 90.46 245.98 60,506 
4B* 9.5 247.12 90.13 247.12 61,068 

(AV=0.30AT) 
When extrapolated to zero pressure: 
V?=61,285+50; V=247.56+0.10 m/sec. 


C.°=20.4+0.2 cal./deg./mole at 90.13°C. 








* The authors believe that this is the lowest pressure for which the 
velocity of sound has been experimentally determined. 


exists between the square of the velocity and the 
pressure, so that extrapolation to zero pressure 
yields a value of the velocity from which the heat 
capacity at zero pressure can be calculated using 
the ideal gas law. 

It will be noticed from the graph that the 
points at the higher temperature are not con- 
sistent for the highest pressures. This is believed 
to be caused by a small degree of polymerization: 
of the butadiene. The determinations were made 
in the order of increasing pressure, and a portion 
of the butadiene remained in the chamber (at 
90°C) for considerable periods of time (different 
in each series) before these high pressure points 
were measured. It was noticed on distilling off 
the material used for series 3 from an auxiliary 
storage bulb at the conclusion of the work that 
a trace of slightly viscous liquid remained. How- 
ever, this affects the lower pressure points only 
slightly, if at all, because except for series 2 (the 
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TABLE III. Heat capacity of gaseous butadiene-1,3 
from velocity of ultrasound and equation of state at 
25.82°C (298.98°K). 











Cocal. Cy cal’ 

Pres- V at deg.~! deg. 

sure 25.82°C mole mole 

cm of in at at 

Point Hg m/sec. Q P 25.82°C 25.82°C 
14 372 223.84 1.026750 0.975735 17.82 17.77 
1B 226 224.85 1.016130 0.985323 17.56 17.53 
1C 630 222.21 1.045913 0.958587 17.94 17.85 
1D 1055 219.53 1.078650 0.929731 18.02 17.87 








TABLE IV. Heat capacity of gaseous butadiene-1,3 from 
velocity of ultrasound and equation of state at 90.13°C 
(363.29°K). 








Cc, cal. 
deg.~! 
mole~! 


at 
90.13° 


20.44 
21.19 
23.57 
20.87 
21.45 
21.15 
21.35 
20.58 


Ce cal. 
deg. 
mole 


at 
90.13°C 


20.46 
21.22 
23.61 
20.88 
21.48 
21.21 
21.37 
20.58 


Pres- V at 
sure 90.13°C 
cm of in 
Point Hg m/sec. Q a 


2A 235 1.0094934 0.991716 
2B 361 1.0146402 0.987257 
2C 590 1.023757 0.979104 
3A 155 1.0062463 0.994542 
3B 398 245.57 1.0161590 0.985940 
3C 760 244.40 1.0312051 0.973023 
4A 301 245.98 1.0121843 0.989380 
4B 95 247.12 1.0038214 0.996657 





246.69 
245.84 
243.90 
246.75 








first at elevated temperature) fresh butadiene 
was used in each series of determinations. 


CALCULATION OF THE HEAT CAPACITY 
BY EQUATION OF STATE 


The following equation of state was obtained 
for butadiene-1,3 by Meyers, Scott, Brickwedde, 
and Rands.! 


1 1,040,000 9.6 
ad |-——— 52} +—, (6) 


y 


RT v 
This equation was used to calculate P and Q in 
Eq. (1). 

The results of these calculations are sum- 
marized in Tables III and IV. In these tables, 
columns 1 and 2 are as before; column 3 is the 
velocity of sound calculated for the common 
temperature of the series; column 4 is Q, as 
defined by Eq. (4); column 5 is P, as defined by 
Eq. (3); column 6 gives the value of C, calculated 
by Eq. (1); and column 7 lists C,° calculated 
from C, and Eq. (6) by the equation 

| (7) 


TT /2°8 1/076 
C.° —_ c+|(—=) +-(— + 
vL\dT?/, 2v\dT?/, 


HARDY, 


AND SZASZ 


The values of C, calculated from the equation 
of state are plotted against pressure in Fig. 2. 
(Point 2C has been omitted from this plot.) 
Extrapolation yields 


C.° at 298.98°K =17.17+0.20 cal. deg.—! mole 
and 
C,° at 363.39°K = 20.26+0.20 cal. deg. mole“. 
These figures should be compared with 
17.33+0.10 cal. deg.—! mole at 298.98°K 
and 
20.4+0.2 cal. deg. mole at 363.29°K 


obtained by using the velocities at zero pressure 
along with the ideal gas law. The agreement of 
the values of C,° calculated from the individual 
points using Eq. (7) along with Eq. (6) with the 
extrapolated values is indicative of the validity 
of the equation of state. There is, however, a 
noticeable trend with pressure. 


CALCULATION OF THE HEAT CAPACITY FROM 
SPECTROSCOPIC AND MOLECULAR DATA 


The heat capacity of gaseous butadiene has 
been calculated using the frequency assignment 
of Aston, Szasz, and Brickwedde.* 

Point group: Cox. 
Carbon skeleton frequencies (in cm~) 


Ag: wi, 1635; we, 1200; w3, 513, 
By: @4, 1590; @®5, 520. 


TABLE V. Calculated and experimental heat capacities of 
gaseous butadiene-1,3 at 298.98°K and at 363.29°K. 











Cp® obs. (1) Cp® obs. (II) Cp obs. (IIT) Cp® calc. 

cal. deg.~! cal. deg.~! cal. deg.) cal. deg.~! 
T°K mole~! mole“ mole~! mole™ 
298.98 19.32 19.16 19.75 18.40 
363.29 22.4 22.45 23.00 21.56 








y (C—H) rocking frequencies 


A,=665, 912, 1010, 
B,=340, 899, 914. 


3J. G. Aston, G. Szasz, and F. G. Brickwedde, ‘‘The 
heat content and free energy of butadiene-1,3 and isoprene 
above 25° C. Equilibria in the dehydrogenation of n- 
butane and n-butenes to butadiene-1,3,”" submitted to the 
Office of the Rubber Director (July 23, 1943). 
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6 (C—H) bending frequencies 


A,=1276, 1303, 1437, 
B,=1285, 1385, 1470. 


v (C—H) stretching frequencies 


A,=3000, 3088, 3000, 
B,,= 3000 (3). 
C—C torsional = 120. 


The calculation was modified by assuming the 
existence of two forms of butadiene—cis and 
trans. Taking AS=0, AH=AF= 1500 cal./mole, 
equilibrium constants for the conversion of cis 
to trans and thence correction terms to the molal 
heat capacities for the fraction converted per 
degree temperature rise were calculated by 
standard equations. 

The experimental and calculated values are 
compared in Table V. In this table, column 1 
gives the temperature; column 2 gives the value 
of C,° obtained by plotting the squares of the 
velocities and calculating C,° by use of these 
limiting velocity squares, while column 3 gives 
the value obtained by plotting the values of C, 
(calculated by the use of the equation of state 
previously discussed) and extrapolating to zero 
pressure. Column 4 gives the value of C,° ob- 
tained by averaging the C,° values iisted in 
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Fic. 2. The heat capacity of butadiene-1,3. @ and left 
scale—t = 25.82°C. O and right scale—t = 90.13°C. 


Tables III and IV, while column 5 gives C,° 
calculated from spectroscopic data. Evidently 
the agreement leaves much to be desired, and 
indicates that the frequency assignment and/or 
the molecular model are not entirely correct. 
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Thermal Isomerization of Butene-1 to Butene-2 


W. WaLTER McCartHy AND JOHN TURKEVICH 
Frick Chemical Laboratory, Princeton University, Princeton, New Jersey 


(Received September 7, 1944) 


The thermal isomerization of butene-1 was investigated in the temperature range of 450- 
550°C using infra-red absorption as a method of analysis. The vibrational spectra of butene-1, 


butene-2, and isobutene are reported. 


OST of the work on the thermal treatment 

of the butenes was primarily concerned 

with their decomposition products,’ polymeriza- 
tion products,? or butadiene synthesis. However, 


1G. Egloff, Reactions of Pure Hydrocarbons, Am. Chem. 
- Monograph Series (Reinhold Publishing Corporation, 
New York, 1937), pp. 333-4. 
? Wheeler and Wood, J. Chem. Soc. 1930, p. 1819. 
’Federov, Smirnova, and Semenov, J. App. Chem. 
(U.S.S.R.) 7, 1166 (1935). 





Runge and Mueller-Conradi* passed butene-1 
through an empty quartz tube at 550°C to 
produce 87 percent of butene-2, no isobutylene, 
and very little liquid polymer. Hurd and Gold- 
sky® studied the pyrolysis of butene-1 and 
butene-2 in the temperature range of 500—700°C 
and contact times of 9-13 seconds. No shift of 


*U.S. Patent 1,914,674, June 20, 1933. 
5 J. Am. Chem. Soc. 56, 1812-15 (1934). 
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the double bond was observed at the lower tem- 
perature. Only at 650°C was there a shift of the 
double bond and this was accompanied by ex- 
tensive decomposition. At 700°C the decom- 
position was so great that the isomerization was 
reduced to several percent. Finally Krause, 
Nemtzov, and Soskina® reported a small amount 
of butene-1 in the products of polymerization of 
butene-2. 


6 Krause, Nemtzov, and Soskina, J. Gen. Chem. 
(U.S.S.R.), pp. 343-55 (1935). 
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A primary difficulty in the previous studies of 
the migration of the double bond is the difficulty 
of analysis for butene-1 and 2. The use of infra- 
red absorption, however, makes this analysis 
very simple. It was, therefore, thought worth 
while to make a study of the thermal isomeriza- 
tion of butene-1 to butene-2. 


MATERIALS AND APPARATUS 


Tank butene-1 was obtained from the Mathe- 
son Company, East Rutherford, New Jersey. It 
was purified by drying over sodium hydroxide 
asbestos. Vapor pressure measurements by Dr. 
R. K. Smith indicated that it was pure butene-1. 

The static system consisted of a reaction 
chamber of Pyrex glass of 200-cc capacity, a 
mercury manometer, a gas burette, and an all- 
glass sampling system. After evacuation of the 
system a volume of butene-1 was measured in 
the gas burette and then transferred to the 
reaction chamber. Changes in pressure during the 
course of the run were used to indicate either 
polymerization or cracking. After a suitable 
length of time, the sample was drawn off and 
condensed in a liquid air trap in preparation for 
analysis. 

The flow system consisted of a reactor which 
was either a Pyrex tube of 8-mm I.D. with a 
20-cm heated length, the same tube packed with 
Fenske glass spirals (3%"), or an iron tube of 
9-mm I.D. with a 20-cm heated length; flow 
meters before and after the reactor, a gas drying 
and a gas sampling system. Temperatures were 
measured by a Chromel-Alumel thermocouple. 
The iron pipe was joined to the Pyrex tubing by 
Picein wax joints. The latter were water-cooled 
to prevent melting. Contact time was calculated 
from t=DV,T2/V27Ti, where D is the time in 
which volume V2 of gas measured at room tem- 
perature 7, flows through the reaction tube of 
volume V; at temperature 7).’ 


ANALYSIS 


The infra-red spectrum of the three butenes is 
given in Fig. 1 and Table I. The spectrum was 
studied from 2-16 mu in a cell 20 cm long and 
at a pressure of 100 mm. The 5.45-mu band of 
butene-1 was used as a measure of the amount of 


7K. Jellinek, Lehrbuch der physikalischen Chemie (Fer- 
dinand Enke, Stuttgart, 1937), Vol. 5, p. 83. 
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TABLE I. 
Butene-1 Butene-2 Isobutene 
Intensity of Intensity of Intensity of 
cm~! absorption cm~! absorption cm~! absorption 
2.31 4329 15 2.30 4348 25 | 4329 25 
3.35 2985 75 3.33 3030 75 3.01 3322 
5.46 1832 35 4.03 2481 10 3.37 2967 75 
6.08 1645 65 4,99 2004 20 4.48 2232 10 
6.82 1466 70 6.03 1658 30 5.59 1789 35 
7.59 1318 40 6.45 1550 6.01 1664 75 
8.53 1172 20 6.88 1453 6.80 1471 85 
9.20 1087 40 7.08 1412 75 7.17 1395 85 
9.98 1010 90 8.81 1135 20 7.73 1294 50 
10.85 922 100 9.64 1037 45 9.48 1055 55 
12.43 805 25 10.15 985 75 10.01 999 35 
15.80 633 40 11.25 889 20 10.96 912 95 
12.10 826 10 11.15 897 95 
14.93 670 100 11.34 882 95 
14.84 674 10 
butene-1 present. Mixtures of known com- Taste II. 
position of the two butenes were made up and ‘ 
a working curve obtained by plotting butene-1 Contact conversion 
P ° ° Run Tube time + ee to butene-1 
against the absorption. The infra-red rocksalt = 
spectrometer employed in this investigation was 4 i empty 32 _ 560 pon 
that described by Barnes, Brattain, and Seitz. B-33 Pyrex packed 4.6 sec. 545 none 
B-34 Pyrex packed 10.0 sec. 550 none 
B-35 iron 8.7 550 22.5 
RESULTS B-36 iron 7.4 550 4.0 
B-38 iron 14.0 550 none 


The results obtained are presented in Table II. 
The run, B-18, at 450°C was carried out in a 
static system at one atmosphere pressure. During 
the first half-hour the pressure decreased slightly 
followed by a subsequent increase of up to 818 
mm at the end of the two-hour run. There was 
substantial decomposition as shown by the brown 
liquid collecting at the outlet tube. An attempt to 
obtain runs with the static system at 550°C 
resulted in such extensive decomposition that an 
adequate sample could not be obtained for 
analysis. 

In the flow system, isomerization was obtained 
with the iron tube and with that, only when it 
was clean. A value of 22.5 percent conversion was 


8 Barnes, Brattain, and Seitz, Phys. Rev. 48, 582 (1935). 








obtained with the clean tube. Subsequent runs 
with the same tube gave increasingly lower con- 
version values. A small amount of carbonaceous 
material was scrapped from the interior of the 
iron tube after these runs. Otherwise, there were 
no other indications of decomposition. 

These results indicate that there is no thermal 
isomerization of butene-1 up to 550°C unless 
there be some concommitant decomposition. 

We wish to express our thanks to Professor 
H. D. Smyth and L. G. Smith for the courtesies 
of the Palmer Physical Laboratory, to the M. W. 
Kellogg Company for financial sponsorhip of 
this investigation, and to Professor Hugh S. 
Taylor for his inspiring interest. 
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Velocity of Sound and Molecular Association in Liquids 


R. T. LAGEMANN 
Department of Physics, Emory University, Emory University, Georgia 
(Received November 5, 1944) 


Two new, empirical criteria for determining molecular association in liquids are pointed out. 
One is the quantity PcR/T¢ in. which R is a function of the velocity of sound; the other is the 
ratio of the velocity of sound in a liquid to the velocity of sound in the gas of the same com- 
pound. In both cases low numerical values of the ratios are qualitative indicators of association. 





LTHOUGH a great many physical proper- 
ties have been used as evidence for 
molecular association in liquids, the velocity of 
sound does not appear to have been studied in this 
connection. It is perhaps of interest then to point 
out two quantities involving the velocity of sound 
which seem to afford a qualitative test for associa- 
tion. These are empirical in nature, and, as is the 
case with other properties, association is ascribed 
to those substances which behave abnormally. 
The first function proposed is 
PcR/Tc=K, (1) 


where Pe is the critical pressure, T¢ the critical 
temperature, and K a quantity varying from 
liquid to liquid, but remainly roughly constant. 
The quantity R, analogous to the molecular re- 
fraction and the parachor, is quite constant for 
any one liquid under changing conditions. Its 
constancy was first pointed out by Rao! who 
defined it as 
R=MV'/d, 

where V is the velocity of sound in the liquid. 
For the calculation of K, values of the critical 
temperature and pressure were secured from the 
International Critical Tables. Most of the values 
of R were obtained from Rao.'? R for water was 
calculated from the data of Bergmann,* and R 
for acetic acid was calculated from the data of 
Smith and Ewing.‘ 

Table I gives values of K calculated by Eq. 
(1) for several substances. It is seen that K is 
roughly constant for most of the normal liquids 
examined. However, it appears to assume de- 
creased values for those substances which are 
commonly believed to be associated. Unfortu- 


1M. Rama Rao, J. Chem. Phys. 9, 682 (1941). 

2M. Rama Rao, Ind. J. Phys. 14, 109 (1940). 

’L. Bergmann, Ultrasonics (G. Bell and Sons, London, 
1938). 

4A. W. Smith and L. M. Ewing, J. Chem. Phys. 7, 632 
(1939). 
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nately data on fused metals and salts and other 
associating compounds are not available to test 
this observation further. Although K contains 
the value of the molecular weight, and it is pos- 
sible to calculate the degree of association, it 
does not seem worth while to do so here. 
Another apparent indicator of association in 
the liquid state is the ratio of the velocity of 
sound in a liquid to the velocity of sound in the 


TABLE I. Calculated values of PcR/Tc. 








Cc Pe 
kK) (Atmos.) R 





Compound ( PcR/Tc 
Acetic acid 594.7 S4.2 629 61 
Methyl alcohol 513.1 78.7 421 65 
Water 647 217.7 207 69 
Acetone 508.1 47 781 72 
n-propyl acetate 549.3 32.9 1211 73 
n-propyl alcohol 536.8 49.9 806 75 


Ethyl alcohol 516.2 63.1 624 76 


n-butyl alcohol 560 48.4 .1004 87 
Ethyl ether 466.9 35.5 1040 79 
Chlorobenzene 632 44.6 1106 78 
Aniline 699 52.4 1078 81 
Carbon tetrachloride 556.2 45.0 944 76 
n-pentane 470.3 33.0 1160 81 
n-hexane 507.9 29.5 1356 79 
n-heptane 539.9 26.8 1545 76 
n-octane 569 24.6 1746 75 
Chloroform 536.0 53.8 803 81 
Benzene 561.6 47.7 979 83 
Toluene 593.7 41.6 1170 82 
Methyl acetate 506.8 46.3 851 78 
Ethyl acetate 523.2 37.8 1037 75 








gas of the same compound. For the determination 
of this ratio data were taken from the literature. 
Because of the marked changes in velocity of 
sound with temperature in both liquids and gases, 
it was decided to reduce all the gas velocities to 
their values at corresponding temperatures of 
0.700 of the critical temperatures, and use ve- 
locities in liquids measured at about 23°C. To 
bring the measured velocities in gases to those 
at corresponding temperatures use was made of 
the correction formula 


Vi/V2=(T1/T2)}, (2) 
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TABLE I]. Ratio of sound velocity in liquids to sound velocity in gases, 








Liquid 


Gas 














Temper- Velocity Temper- Velocity ’ Viiquia 
ature of sound Ref- ature of sound Ref- — 
Compound (°C) (m/sec.) erence (°C) (m/sec.) erence Vegas 
Water 24 1494 B 179.8 519.6 y 2.88 
Methy] alcohol 23.8 1130 B 86.1 350.5 D 3.22 
Ethyl alcohol 23.5 1207 B 88.2 274.2 D 4.40 
n-butyl alcohol 23 1315 B 119.0 235.9 D 5.57 
Acetic acid 25 1140 E 143.2 210.9 5.41 
Carbon disulphide 23 1149 B 109.1 216.4 A 5.31 
n-pentane 18 1052 B 56.1 195.5 8 5.38 
Ethyl acetate 23.5 1187 B 93.1 213.1 D 5.57 
Acetone 22:5 1203 B 82.6 213.5 fy 5.64 
Methylene chloride 23.5 1064 B 89.6 188.4 ¢ 5.65 
Chloroform 23.5 1001 B 102.1 171.5 D 5.84 
Carbon tetrachloride 23 929 B 116.2 158.2 D 5.87 
Benzene 23 1310 B 120.0 219.4 D 5.97 
n-hexane 23 1113 B 82.4 184.9 Cc 6.02 
4 International Critical Tables. B Bergmann, Ultrasonics (G. Bell and Sons, London, 1938). 
ys 1.C.T. and correction by Eq. (2). Bergmann and correction by Eq. (2). 
* A. W. Smith and L. M. Ewing, J. Chem. Phys. 7, 632 (1939). 
where temperatures are in degrees Kelvin. That substances. Furthermore, the two proposed 


this would correct the data with adequate pre- 
cision was determined by testing it on data for 
several gases found in the I.C.T. Irons® has also 
shown that it applies to several gases over a wide 
range of temperature. 

Table II gives the velocity of sound in several 
gases corrected to 0.700 of the critical tem- 
peratures. In column 3 are recorded velocities in 
liquids measured at the temperatures given in 
column 2. The last column gives the ratio of the 
velocity of sound in the liquid to that in the gas 
of the same compound. It may be seen that for 
those compounds commonly accepted as being 
associated the ratio is smaller than for the 
normal liquids. The relatively high value for 
acetic acid may be caused by the fact that it is 
more associated in the gaseous state than are the 
other compounds listed. According to Laplace’s 
equation, the increased molecular weight would 
lower the velocity in the gas below its value for 
the monomer and raise the ratio above its ex- 
pected value. This explanation is supported by 
the case of oxygen where it has been found that 
the presence of ozone decreased the velocity of 
sound in the gas. 


COMPARISON WITH OTHER METHODS 


It is difficult to compare these results with 


previous methods for determining association 
because the earlier methods do not agree even as 
to the relative order in which to place various 


°E. J. Irons, Phil. Mag. 3, 1274 (1927). 


schemes do not possess high correlation. 

It should be pointed out that the foregoing 
functions of sound velocity cannot be used to 
predict correctly changes in the association of a 
liquid with changes in temperature. In fact, the 
ratio of the liquid to the gas velocity for any 
substance except water (where velocity in the 
liquid increases with increase in temperature) 
indicates an increase in association with increase 
in temperature, in direct contradiction to ac- 
cepted belief. Therefore, these two quantities are 
considered to be only qualitative criteria for 
association. 

It is usually found that increase in the parachor 
with increase in temperature is indicative of as- 
sociation. The same does not appear to be true 
for the analogous function R. From the work of 
Rao? R may be seen to increase slightly with 
increase in temperature in certain non-associated 
as well as certain associated liquids. 

An advantage of the sound velocity criteria for 
association over some other properties is that the 
former are independent of measurements con- 
nected with the liquid surface where unusual 
conditions exist. On the other hand, the well- 
known use of ultrasonics both as a dispersive and 
as a coagulating agent may change the degree of 
association even as the velocitiés are being 
measured. As more data become available it 
would be of interest to test the above considera- 
tions on other substances, particularly some of 
the lower aldehydes and nitriles. 
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Letters to the Editor 








ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 














and 


B = In potA /v¢?, 


bears little formal resemblance to the Brunauer-Emmett- 
Teller (BET) isotherm, 


- Combo (4) 


—p)(pbo—p+cp)' 





. 


third of the month. Because of the late closing date for the sec- ) 
tion no proof can be shown to authors. The Board of Editors vet the two equations agree over a wide range of conditions. | 
does not hold itself responsible Jor the opinions expressed by 4 preliminary step toward explaining this relationship has 
the correspondents. Communications should not in general ex- yeen made. The author! has shown that at least in one 
ceed 600 words in length. case Eq. (4) combined with Eq. (2) gives a result in the 
form =b—age, for a considerable pressure range. In view 
of the recent statement by Harkins and Jura’ that the 
agreement between the two isotherms “‘is extremely re- ) 
The Relationship Between the Brunauer- markable when the very great differences between the two é 
Emmett-Teller Adsorption Isotherm Ho one erent and ae ont yeas a are 
taken into consideration,”’ it shou of interest to pursue 
and the New Isotherm of Jura this investigation further. 
and Harkins Although it has been said that the Jura-Harkins iso- 
H. K. Livincston therm is thermodynamic in origin, this is only because it 
Jackson rn 4g Faoypai and Company, involves the Gibbs equation (2). The isotherm owes its 
Sunteniber 28, 1944 success to the experimental fact that the surface pressure 
(w) varies as a straight-line function of the molecular area 
HE first adsorption isotherm to give good agreement (¢), Over a considerable pressure range and for a large 
with experimental data from all types of vapor-on- umber of systems. If the BET isotherm is combined with 
solid adsorption was that of Brunauer, Emmett, and  §a. (2), the result is: 
a= (RTvm/VZ) In [((po—p+cp)/(po— p) J. (5) 
TasLe I. ) 
Furthermore, . 
Relative pressures at which the 
Constants of straight- equation fits the data o= ZV (po— p) (po— p+cp)/Novmcppo. (6) 
line Eq. (7) within 3 percent 
: wad ” DEES en See See Therefore, from Eqs. (5) and (6) it should be possible to 
2 0.165 0.545 0.20 0.40 determine the conditions where the linear relation between 
= = HH rr ~ x and o applies. This has been done by choosing values of 
25 1.28 2.07 14 _ 60 p/po throughout the pressure range and also various values 
a = — = = of c, and solving Eqs. (5) and (6). It was then possible to 
250 2.30 3.62 07 50 determine the range where Eq. (7) applies. ? 











a 
(VE/RTvm)x = b’ — (Novm/V=Z)me, (7) ‘ 
: e = b V=/R Tim >, n= A /20m?. ae | c 
Teller.’ Jura and Harkins? have developed a new isotherm of 
which gives as good agreement with experimental data 6’ and m are the intercept and slope of the straight line, od t 
up to b/po=}, and they find even better agreement at respectively.) The results (Table I) show that Eq. (7) a 
higher pressures. Their isotherm was obtained by dif- [Lor Eq. (1)] is in good agreement with the BET isotherm ; tl 
ferentiating the pressure-area relation for condensed over a considerable pressure range. The range is especially ) 
monolayers: broad if c is from 10 to 250. Experience to date indicates b 
that most systems have c values in this range.* 5 Naturally, “4 si 
™=b—ao (1) plots of In p vs. vm?/v?, using the values of v obtained from i 
and combining with the Gibbs adsorption isotherm for 4: (4), also give straight lines over the pressure ranges 2 
solid surfaces [Eq. (678) of the ‘Equilibrium of hetero- indicated by Table I. ; $ 
geneous substances], which, using modern symbols, is: We may conclude that the observed correlation between ie 
the Jura-Harkins and BET isotherms is in agreement with in 
dx =(RT/VZ)(v/p)dp. (2) the results obtained from a mathematical treatment of the ec 
(Note: The synibols used are identical with those of Har- two isotherms. 4 
kins and Jura.’) The resultant equation 1 Brunauer, Emmett, and Teller, J. Am, Chem, Soc. 60, 309 (1938); * 
Brunauer, The Adsorption of Gases and Vapors. Vol. I, Physical Ad- a di 
In p= B-—A/v? sorption (Princeton University Press, Princeton, 1943). a d 
’ 2G. Jura and W. D. Harkins, J. Chem. Phys. 11, 430 (1943). x. e 
where tf e ae and by Jura, J. -. 7. Ly" 6! 1366 (1944). | : et 
ivingston, m m : 
A =a2*V?/2NoRT (3) ’ Gans, Brooks, and Boyd, Ind. Eng. Chem. Anal Ed. i4, 396 (1942). tri 
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A Relation Between Bond Multiplicity and 
Interatomic Distance 


J. L. KAVANAU 
California Institute of Technology, Pasadena, California 
November 2, 1944 


N this paper an empirical equation relating the inter- 
atomic distances of covalent bonds to the bond mul- 
tiplicity is advanced. 

Pauling, Brockway, and Beach, and Penney,? have 
treated aspects of this problem using known interatomic 
distance values for carbon single, double, and triple bonds 
and for the bonds in graphite and benzene. The former 
investigators assumed bond-linkage values of 4/3 and 3/2 
for graphite and benzene bonds, while Penney by a quan- 
tum-mechanical calculation determined bond-order values 
of 1.45 and 1.623 for these bonds. 

Certain simple considerations suggest an equation of the 
following form, involving the bond multiplicity N for 
single, double, and triple bonds, 


D=a+bC%. (1) 


D is the interatomic distance in angstrom units, a and b 
are constants in angstrom units determined by the specific 
atom pair, and C is a dimensionless constant. The value of 
the constant C determined by applying this equation 
suggests that it is equal to the function (m:+m2—1)/ 
(nm: +m2.+1), where m, and m2 are the principal quantum 
numbers of the valence electrons in the bonded atom pair. 
Thus, I have modified Eq. (1) to 


ni+ne— 1\% 
- eee 6 2 
- a Coe (2) 
By use of the equation 
_ 
log (D—a) = N log (ate) +1og b, (3) 


a linear semilog plot of (D—a) as a function of N may be 
made. Since (m:+n2.—1)/(mi+m2+1) is a dimensionless 
constant which determines the slope of the semilog plot, 
the curves will be of the same slope for all bonds between 
atoms homomorphic in the principal quantum number of 
the valence electrons. 

For convenience in treating single, double, and triple 
bonds, the parameters a and 6 may be eliminated by 
simultaneous solution of Eq. (2) in N=1, 2, 3 giving 


D,—D3_m+m—1 
D,—-D, m+m+l1 





(4) 


The principal quantum number for the valence electrons 
in carbon and nitrogen is 2, so that (m1-++-m2—1)/(m1+m+1) 
equals 3/5 for such atom pairs. 

For carbon bonds the single-bond distance is 1.54 
+0.02A,3 and is known to a high degree of accuracy in 
diamond, being 1.542A.4 For the double bond the various 
determinations of 1.325+0.005A® and 1.353+0.01A® in 
ethylene and 1.330+0.005A° in allene are known. For the 
triple bond the value for acetylene is 1.204+0.005A.’ 
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Using the single- and triple-bond distances, since there is 
some doubt about the double-bond value, Eq. (4) yields 
a value of D,=1.330A. This checks well with the experi- 
mental values for the double-bond distance. The values 
for the parameters are a= 1.015 and }=0.875. 

For nitrogen bonds the value for the single bond in 
hydrazine is 1.47+0.02A.8 The double-bond value for 
transazobenzene is 1.23+0.03A,9 and the value for the 
triple bond in the nitrogen molecule is 1.094+0.005A.’ 
Using the more accurate single- and triple-bond values, 
Eq. (4) yields D2=1.235A which checks closely with the 
experimental value. From Eq. (2) the values for the 
parameters are a=0.883 and b=0.979. 

Using the parameters determined for carbon bonds, Eq. 
(2) becomes 


D=1.015+0.875(3/5)%. (5) 


This equation gives the graphite and benzene interatomic 
distance values of 1.42A and 1.39A if N assumes the 
values 1.508 and 1.659, respectively. These check fairly 
closely with Penney’s values of 1.45 and 1.623 calculated 
for the bond order. Accordingly, it appears that NV in my 
equation has the same significance for bonds of fractional 











TABLE I. 
De 
Atom Di Ds Covalent D2 
pair Exp. Exp. radii Eq. (4) 
C-N 1.47 +0.02A¢ 1.154 +0.005A° 1.265A 1.272A 
P-P 2.21+0.02A¢ 1.86 +0.005A¢ 2.00A 2.006A 








(1936) O. Brockway and H. O. Jenkins, J. Am. Chem. Soc. 58, 2036 
®See reference 7 of text. 
¢L. R. Maxwell, S. B. Hendricks, and V. M. Mosley, J. Chem. Phys. 
3, 698 (1935). 
4G. Herzberg, Ann. d. Physik 15, 677 (1932). 


multiplicity as Penney’s bond order. The values 1.508 and 
1.659 are only eight-thousandths removed from the frac- 
tions 3/2 and 5/3, which are each 1/6 greater than the 
respective linkage values assigned by Pauling, Brockway, 
and Beach, this being a contribution of the resonance 
energy. 

I know of no other atom pairs for which the single-, 
double-, and triple-bond distances are known which would 
permit further checking of Eq. (2). The single- and triple- 
bond distances are, however, known for carbon-nitrogen 
bonds and phosphorus bonds so it is interesting to predict 
the respective double-bond distances and compare them 
with the distances indicated in Pauling’s table of covalent 
radii.” The results are tabulated in Table I. 
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